arXiv:math-ph/0503032v2 23 Mar 2005 


Quantum Chaos: Spectral 
Analysis of Floquet Operators 




THE UNIVERSITY OF 

MELBOURNE 


Quantum Chaos: Spectral 
Analysis of Floquet Operators 


James Matthew McCaw 


Submitted in total fulfilment of the requirements 
of the degree of Doetor of Philosophy. 


School OF Physics • The University of Melbourne 
December 2004 Australia 



Copyright © James Matthew MeCaw 2004-2005. 

This work is licensed under the Creative Commons Attribution-NoDerivs 2.0 
License (Australia) http://creativecommons.Org/licenses/by-nd/2.0/au/ 

You are free to copy, distribute and display this work, and to make commercial use 
of this work, provided that you give credit to the author James McCaw. You may 
not alter, transform, or build upon this work. For any reuse or distribution, you 
must make clear to others the license terms of this work. Any of these conditions 
can be waived if you get permission from the copyright holder, James McCaw. 
Your fair use and other rights are in no way affected by the above. 


Online edition. 


National Library of Australia Cataloguing-in-Publication entry: 
McCaw, James Matthew. 

Quantum Chaos: Spectral Analysis of Floquet Operators. 

Bibliography. 

ISBN 0 7340 3071 1. 

ISBN 0 7340 3068 1. (hardback print edition) 

1. Quantum chaos. 2. Differential operators. 3. Floquet theory. 

I. University of Melbourne. School of Physics. II. Title. 

530.12 


The online version of this thesis is archived at the University of Melbourne ePrints 
Repository and is part of the Australian Digital Theses project. It has the Open 
Archives Initiative identifier oai:unimelb.edu.au:848, and is accessible at 
http://errol.oclc.Org/oai:unimelb.edu.au:848.resource 


To Victoria. 



Abstract 


The Floquet operator, defined as the time-evolution operator over one period, 
plays a eentral role in the work presented in this thesis on periodieally perturbed 
quantum systems. Knowledge of the speetral nature of the Floquet operator gives 
us information on the dynamies of sueh systems. The work presented here on 
the speetrum of the Floquet operator gives further insight into the nature of ehaos 
in quantum meehanies. After diseussing the links between the speetrum, dynam¬ 
ies and chaos and pointing out an ambiguity in the physics literature, I present a 
number of new mathematical results on the existence of different types of spectra 
of the Floquet operator. I characterise the conditions for which the spectrum re¬ 
mains pure point and then, on relaxing these conditions, show the emergence of a 
continuous spectral component. The nature of the continuous spectrum is further 
analysed, and shown to be singularly continuous. Thus, the dynamics of these sys¬ 
tems are a candidate for classification as chaotic. A conjecture on the emergence 
of a continuous spectral component is linked to a long standing number-theoretic 
conjecture on the estimation of finite exponential sums. 
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Chapter 1 


Introduction 


Contents 


1.1 Summary of results. 3 


The main work presented in this thesis is an extension of a number of math- 
ematieal results concerning the spectral analysis of time-evolution operators in 
quantum mechanics for a class of Hamiltonian systems. While the results are of a 
technical nature and do not lend themselves to a simple “physical” interpretation, 
they do have potential impacts on the broad field of quantum chaos; specifically 
the link between the spectrum of an operator^ and the dynamics of the system. 

While this link to quantum chaos turns out to have little impact on the core 
results presented, it sets the context for the research and drives the course of work 

'The operator of interest is either the Hamiltonian or the time-evolution operator over one 
period, depending upon the type of system under consideration. 
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Introduction 


presented in the two main ehapters. I will begin with an investigation of the so¬ 
lutions to the Sehrodinger equation when the Hamiltonian is some simple system 
(e.g., the harmonic oscillator) which is perturbed in a periodic fashion. I will then 
characterise the spectrum and extend a number of known results in the field of 
spectral analysis. Along the way, I will build unitary equivalents of a number 
of well known self-adjoint theorems from functional analysis and discover links 
between physical conjectures in the physics literature and deep number-theoretic 
conjectures. The results presented in Chapter 6 and Chapter 7 constitute the main 
body of research undertaken. 

Chapter 2 provides an introduction to both classical and quantum chaos, out¬ 
lining definitions, techniques and model systems widely used in the literature. 
Chapter 3 introduces the Floquet operator, an essential tool in the investigation of 
time-periodic quantum mechanical systems. 

To discuss the links between dynamics and spectra presented in Chapter 5 
properly, an understanding of functional analysis and the mathematics of Hilbert 
spaces is required. Chapter 4 provides this background, and pointers into the 
mathematics literature where required. This chapter also develops the background 
theory and notation used in Chapter 6 and Chapter 7. 

As already mentioned. Chapter 5 constitutes a more detailed examination of 
the literature that motivated my study of the spectrum of operators. I discuss the 
link between quantum dynamics and spectral analysis of the Hamiltonian oper¬ 
ator (time-independent systems) or the time-evolution operator (time-dependent 
systems). I will also identify a general ambiguity in the literature, pointing out the 
pitfalls and also a potential resolution to this problem. 



1.1 Summary of results 


3 


1.1 Summary of results 

The results of this thesis are presented in detail in Chapter 6 and Chapter 7. I 
eonsider systems deseribed by the time-dependent Hamiltonian 

OO 

H{t) = Ho + \wJ2 

n=0 

where Hq is the time-independent Hamiltonian for some simple system (e.g., the 
harmonic oscillator), W is an operator describing the perturbation and A is a 
strength parameter. The spectral properties of the time-evolution operator over 
one kick period, T, are investigated. 

In Chapter 6 I show that, for arbitrarily strong perturbation strengths A, the 
spectrum of the time-evolution operator remains pure point if the perturbation 
operator W satisfies the condition 

< OO. 
n 

The (pn are the basis states of the Hq system. Essentially, we need when 

written in terms of the basis states of Hq, to be li convergent.^ This is a non- 
perturbative result and indicates that the system described by H (t) is stable—the 
dynamics does not change in a fundamental way due to the perturbation. 

The condition on W is relaxed in Chapter 7 and the possibility for fundamen¬ 
tally different dynamical behaviour is shown to arise in the case where Hq is the 
harmonic oscillator, with eigenvalues (ignoring the l/2hiu term) given by nhuj. 
When the ratio between the kicking period, T, and the natural frequency, uj, of the 
unperturbed system is irrational the Floquet operator is shown to obtain a continu¬ 
ous spectral component. This result is proved for rank-iV perturbations (i.e., W is 

^To be a well defined state in the Hilbert space, it must be I2 convergent. That it is also li 
convergent is a further restriction which a state in the Hilbert space may or may not satisfy. 
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a rank-iV operator built from N projections), extending the previous rank-1 work 
of Combescure [28]. 

The analysis is then extended to the case where Hq is not the harmonic oscilla¬ 
tor, but some general pure point system with eigenvalues described by an arbitrary 
order polynomial in n (rather than a first order polynomial as is the case for the 
harmonic oscillator). I show that the question of whether a continuous Floquet 
spectrum will exist is equivalent to a number-theoretic conjecture presented by 
Vinogradov [112] over fifty years ago. The conjecture concerns the estimation of 
Weyl sums—finite sums of the exponential of arbitrary order polynomials. The 
optimal estimation of such sums for the case where the coefficients of the polyno¬ 
mial are rational is known. For the case where the coefficients are irrational, the 
current best estimates are not believed to be optimal. Vinogradov’s conjecture is 
that the optimal estimation for such sums is that of the rational coefficient case. 

If this conjecture were shown to be true, that is, if such Weyl sums were shown 
to have a bound given by the rational coefficient case, then the strongest condi¬ 
tions for the existence of a continuous spectrum for the Floquet operator could be 
established. 

The final piece of work presented in Chapter 7 concerns the further classifi¬ 
cation of the continuous spectrum of the Floquet operator. Extending the results 
of Milek and Seba [90], I show that the continuous Floquet spectrum is in fact 
singularly continuous. By the arguments presented in Chapter 5, it is clear that 
the existence of a singularly continuous component of the Floquet operator is a 
necessary condition for the time evolution of such systems to show chaotic struc¬ 


ture. 
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Chaos in classical and quantum mechanics 


The classical study of chaos is now a well established and flourishing field of 
research in mathematics and mathematical physics. Chaotic behaviour seems to 
pervade a large number of dynamical systems, and an appreciation of it is essential 
in many contexts, both applied and theoretical. The classic example of the earth’s 
weather patterns always springs to mind when chaos is mentioned. 

The microscopic world, however, is not governed by the laws of classical dy¬ 
namics. In the realm of small quantum numbers the dynamics of a system is 
governed by the Schrodinger equation. In such systems, the simple and elegant 
definitions of chaos such as positive Lyapunov exponent, which hold for classical 
systems, are not applicable. In fact, there is no universally accepted definition for 
what constitutes chaos in quantum mechanics. The realisation that the concept 
of a phase space is not easily translated to quantum mechanics (the uncertainty 
relation makes it impossible to assign both a position x and a momentum p si¬ 
multaneously to a quantum object) immediately precludes a simple link between 
classical characterisations of chaos and quantum characterisations of chaos. The 
route past this difficulty is not obvious, and much of the work so far has been 
centred around finding appropriate tools. As yet, no consensus has been reached 
on an appropriate definition of quantum chaos. 

That said, from the early 1980s onwards many papers have been published that 
adapt classical concepts and bring them into the quantum world. There are also 
papers that introduce completely new concepts not based on equivalent classical 
ideas, in an attempt to make sense of quantum dynamics. I will survey some of 
these ideas in this chapter. 
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Three broad questions motivate the eontinued investigation of quantum ehaos. 

(a) Are there quantum systems that display chaotic behaviour at the quantum 
level? 

(b) What properties of a quantum mechanical system determine whether or not 
the corresponding classical system will display chaotic behaviour? 

(c) Do we come to the same conclusions irrespective of the choice of defini¬ 
tional criteria for quantum chaos? 

The first question is very difficult to answer, in part because there is no clear 
meaning to the expression “chaotic behaviour” when dealing with microscopic 
phenomena. I will have more to say about this later. 

The second question is intimately linked to the Correspondence Principle and 
theories of quantum measurement. It is a generally held belief of physicists (al¬ 
though some do not agree, see [83]) that the classical world is in fact quantum 
mechanical—how we obtain the classical world from quantum mechanics is a 
topic of extensive investigation, and needless to say, this area of fundamental 
physics is infamous for its interpretational difficulties and seemingly inconsis¬ 
tent behaviour. Mixing chaos into the discussion can only make matters more 
interesting. 

The final question is not independent of the previous two, but important to con¬ 
sider when examining the literature. Many attempts to look at particular aspects 
of the dynamics of quantum systems have been made. They are based around cer¬ 
tain toy systems and models. Some papers, as a consequence of their definition 
of quantum chaos, come to the conclusion that there is no such thing [63]. That 
is, they conclude that quantum systems do not display chaotic behaviour. Other 
papers, simply as a result of a different starting point, come to the opposite con¬ 
clusion. One can certainly say that the pursuits of researchers has uncovered rich 
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variations in the dynamics of simple quantum systems, some of which are now 
accessible to experiment. New experiments should provide great insight into the 
numerous models and ideas currently being discussed in the literature. The final 
word on this exciting field of research is a long way off. 

I will now proceed to give an overview of the techniques used in classical 
and quantum chaos research, and discuss model systems that receive the most 
attention from the research community. The discussion leads naturally to the new 
work presented in this thesis on the spectral analysis of unitary operators. 

2.1 Chaos in classical systems 

The theory of chaos in classical mechanics is well established and understood. 
The most accessible definition of chaos is obtained from the phase space approach 
and the growth of the largest Lyapunov exponent for a dynamical system. The path 
drawn out in phase space {x, p) for a particle is determined by solving the equa¬ 
tions of motion. If initially similar phase space points diverge exponentially as 
time progresses, a Lyapunov exponent for the system is positive, and we say that 
the system is chaotic. The system displays sensitivity to initial conditions. Con¬ 
versely, systems which only experience power law type separation of trajectories 
do not have any positive Lyapunov exponents and are deemed non-chaotic. 

The Lyapunov exponent definition of chaos is only one of the available char¬ 
acterisations. The famous paper Period Three Implies Chaos by Tien-Yien Li and 
James A. Yorke [87] shows that if a system described by a continuous map^ has 
a period of three^ then it will have periods of all integers and, importantly, also 
have an uncountable number of initial states which never come even close to being 

'While Li and Yorke’s original paper was concerned with interval maps, it has been extended 

to the general frame of topological dynamics. See [9]. 

^i.e.. After three iterations of the dynamics, we return to the exact initial point 
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periodic. That is, the evolution of those states never returns to a point arbitrarily 
close to the initial state. 

A third way to characterise chaos in classical systems is to quantify the infor¬ 
mation required to describe the dynamics of a system. Essentially, if a dynamical 
system requires an exponentially increasing number of bits to accurately simu¬ 
late its evolution for a linearly increasing time, then the information content is 
high, and we say the system is chaotic. If the increase in bits required follows 
only a power law as the desired time of simulation increases, the system is non- 
chaotic. Information theory approaches to chaos also include considerations of 
entropy production in systems. For a beautiful discussion of the links between 
unpredictability, information and chaos, see [23]. 

The different characterisations of chaos presented here are believed to be es¬ 
sentially equivalent; attempting to establish this equivalence is an active field of 
research. As an example, Blanchard et. al. [9] have recently shown that the posi¬ 
tive entropy characterisation of chaos implies Li-Yorke chaos. For a review of the 
different approaches to chaos see, for example, the review article by Kolyada and 
U. Snoha [81]. 

In classical mechanics, the simplest model systems in which we observe chaos 
are iterative, rather than continuous. The logistic map and the Henon-Heiles sys¬ 
tems display chaotic behaviour for certain parameter values. Of the continuous 
systems, the double pendulum is arguably the most elegant example of a chaotic 
system. The equations of motion are easily solvable, yet the resulting dynamics 
shows the complexities of chaotic behaviour. Importantly, note that the system is 
conservative—the total energy is constant over time. 
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2.2 Kicked systems as a model testing ground for 
chaos 

The field of non-eonservative dynamics is a rich one, and chaotic behaviour is a 
common feature. In particular, it turns out that systems perturbed by an external 
sharp pulse (a “kick”) have interesting dynamical properties. For example, con¬ 
sider a spinning top—a childhood toy. If left to its own devices, it will behave 
quite predictably, its main axis of rotation processing slowly, smoothly and pre¬ 
dictably. There is certainly no chaotic motion. If, however, the spinning top is 
periodically kicked by a short sharp pulse, the motion turns out to be quite unpre¬ 
dictable; the system is chaotic [57]. 

This observation has motivated many researchers in the field of quantum chaos 
to consider quantum systems which are periodically perturbed. They present an 
attractive combination of physical realisability and mathematical tractability. Fea¬ 
turing strongly in the literature are the kicked rotator (work by Izrailev and Shep- 
elyanskii [71], Grempel and Prange [51], Casati et. al. [22], and Dittrich and Gra¬ 
ham [37]) and the kicked harmonic oscillator (see Combescure [28], Graham and 
Hiibner [50], and Daly and Hefernan [31]). The quantum kicked top is also exam¬ 
ined in great detail by Haake et. al. [57,58]. A whole array of other kicked systems 
have been considered in relation to the study of quantum chaos—see [4,5,91,92] 
for just some examples. Many of the papers referred to in the next section are also 
concerned with the analysis of kicked systems. 

2.3 Characterisations of quantum chaos 

As already mentioned, there is no simple way to take the classical characterisa¬ 
tions of chaos and use them in the quantum context. As a consequence, a great 
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deal of research has been done attempting to find appropriate characterisations of 
quantum systems which provide a clear link to the classical concepts of chaos. In 
this section, I will review a number of these attempts, highlighting the diversity in 
the field of quantum chaos. While I separate the discussion into broad sections, 
this is somewhat artificial. There are significant overlaps between some of the 
fields. Many of the papers referred to bring together many aspects of quantum 
chaos. 

The paper by Caves and Schack [23], concerning information-theoretic ap¬ 
proaches to dynamics, highlights a number of the issues that make the study of 
chaos so much more difficult in quantum systems, as opposed to the classical 
systems. 

2.3.1 Gutzwiller, periodic orbits and random matrix theory 

The field of quantum chaos was arguably born with the work of Gutzwiller [55]. 
His semi-classical work on calculating the energy eigenvalues for Helium via clas¬ 
sical periodic orbits identified the importance of the chaotic nature of the classical 
dynamics. For more details on his work, see [54, 56]. An excellent review was 
conducted by Heller and Tomsovic [59]. A few examples of the influence of 
Gutzwiller’s work in the literature is seen in the papers of Eckhardt, Wintgen et. 
ah, Tomsovic and Heller, and Keating [39,78,110,113]. 

The ideas of Gutzwiller are still of central importance. In an attempt to get 
to the essence of his work, the study of quantum billiards has been developed. 
Classically, the motion of a particle in a two-dimensional bounded space (such as 
a billiard table) may be chaotic. The motion is intimately linked to the geometry 
of the box. By considering a quantum particle constrained to such a geometry, 
and using Gutzwiller’s ideas relating the classical periodic orbits to the quantum 
energy levels, one can gain great insight into the dynamics of the quantum system. 
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Importantly, these investigations have also lead to aetual experimental investiga¬ 
tions in microwave cavities, such as those by Sridhar and Lu, presented in [108] 
and references therein. 

These considerations are intimately related to another technique of great im¬ 
portance in the field of quantum chaos— random matrix theory. The eigenvalue 
statistics of non-chaotic and chaotic systems are, in a sense, generic, by which I 
mean that the energy eigenvalues of chaotic systems typically have the same sta¬ 
tistical structure. They are accurately described by random matrix theory—the 
field of mathematics concerned with the statistics of eigenvalues of classes of ran¬ 
dom matrices. Billiard systems are an essential tool in work on random matrix 
theory [13]. See the work of Smilansky [107], Kettemann et. al. [79] and the 
proceedings [26] along with references therein. 

2.3.2 de Broglie Bohm theory 

Since the initial conception of quantum mechanics, there have been persistent at¬ 
tempts to find a “classical” interpretation of the theory. The dominant Copenhagen 
interpretation of quantum mechanics is most certainly not such an interpretation. 
The de Broglie Bohm interpretation of quantum mechanics [62], largely ignored 
in the physics community, attempts to provide a strong conceptual link between 
quantum and classical dynamics; it is based around the Hamilton-Jacobi equa¬ 
tions of motion. A quantum particle is just like a classical particle, except that the 
classical potential it exists in is supplemented by a quantum potential. 

The quantum potential is given through the solution of the Schrodinger equa¬ 
tion, and thus, results in de Broglie Bohm theory are equivalent to those in stan¬ 
dard quantum mechanics.^ Given that a quantum particle now has a trajectory, a 
^For over seventy years, this view of the equivalence has prevailed. This is in fact one of the 
reasons the theory is usually dismissed—it is often seen as nothing more than a rewrite of the 
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phase space picture of the dynamics is realisable. Thus, both a quantum Lyapunov 
exponent and hence quantum chaos are definable just as in classical mechanics. 
The quantum potential blurs the trajectory somewhat, but one can think of a “flux 
tube” in phase space, describing the evolution of the system. 

Papers employing the de Broglie Bohm theory to examine quantum chaos in¬ 
clude those by de Alcantara Bonfim et. al. [33], de Polavieja [35], Schwengelbeck 
and Faisal [42,102], and Wu and Sprung [114]. The de Broglie Bohm theory also 
plays a central role in a number of the other ideas listed in the following sections. 

2.3.3 Open quantum systems 

The vast literature on quantum mechanics in open systems also provides a strong 
link to classical chaos. Tiny environmental interactions which have no effect 
on the classical systems turn out to be strong enough (easily!) to completely 
change the quantum dynamics. The coherent interference effects responsible for 
the “quantum suppression of chaos”'* are destroyed and the quantum systems show 
rather classical behaviour, including chaotic behaviour. Fritz Haake’s book Quan¬ 
tum signatures of Chaos [58] and the many references therein provide an in depth 
review and excellent discussion of these phenomena. The ideas of Zurek, further 
discussed by Paz, have also been influential [93,119-121]. There has also been 
experimental work in this area, including the early work of Bliimel et. al. [ 10,11 ]. 
On the theoretical side, of interest is much of the work by Cohen, Dittrich et. al., 
Grobe et. al, and Kohler et. al. [25, 37, 38, 52, 80]. The links back to random 

normal quantum theory. However, in recent years there have been some, such as Ghose [47], who 
claim there are incompatibilities between the theories. Experimentally realisable tests have even 
been suggested by Golshani and Akhavan [48]. The work is, however, highly controversial and 
disputed. Struyve and De Baere [109] and other authors referenced therein provide arguments 

against this suggested incompatibility. 

will detail the meaning of this expression in Section 2.4. 
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matrix theory and level statistics are strong. This work also provides significant 
insight into Question (b) on page 7. 

2.3.4 Other interesting techniques 

Here I briefly list a number of other ideas which have been put forward in the liter¬ 
ature along with a number of references. Again, many of the ideas draw upon the 
basic ideas of random matrix theory, de Broglie Bohm theory and open quantum 
systems. 

The Loschmidt echo 

A very interesting approach to unifying the concept of chaos in classical and quan¬ 
tum mechanics has been presented by Jalabert and Pastawski [72] and further re¬ 
viewed by Cucchietti et. al. [29]. In classical mechanics, small perturbations in 
the equations of motion lead to the exponential divergence of trajectories in phase 
space. It turns out that in quantum mechanics, while small changes to initial con¬ 
ditions do not lead to significantly different dynamical behaviour, small changes 
in the Hamiltonian can lead to significant variations in the time evolution. The 
overlap of initially identical wave functions is measured and a Lyapunov expo¬ 
nent is extracted. This recent idea appears promising as it provides an opportunity 
to directly compare chaotic structure in classical and quantum systems. 

The Quantum action 

The quantum action was introduced by Jirari et. al. [74,75] and utilised by Caron 
et. al. [17-19]. The aim is to unify the characterisation of chaos in classical and 
quantum mechanics by introducing a quantum action analogous to the classical 
action. Jirari et. al. conjecture that [75] (quote) 
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For a given classical action S with a local interaction V(x) there is a 
renormalized/quantum action 

J dt—x^ — V{x), 

such that the transition amplitude is given by 

C^(^final) ^final) ^ CXp 

where Xdass denotes the classical path corresponding to the action S. 

Note that the mass term, fh, and the potential, V, in the action are quantum param¬ 
eters. The quantum action takes into account quantum corrections to the classical 
motion. Once obtained, the tools of classical mechanics may be applied as the 
mathematical form corresponds exactly with the classical action. The integral is 
taken only over the classical path. 

This direct link allows the definition of chaos in terms of the classical action 
to be taken over to the quantum dynamics directly. 

It must be noted, as acknowledged by Jirari et. al. [74], that there is no proof 
of the conjecture. Numerical evidence is presented by Jirari et. al. and Caron et. 
al. that indicates that it seems to be reasonable in a range of cases. 

Entropy approaches and information theory 

As stated by Slomczyhski and Zyczkowski in [106], “the approach linking chaos 
with the unpredictability of the measurement outcomes is the right one in the 
quantum case”. To measure this unpredictability, they introduce a generalised 
notion of entropy. The inclusion of the measurement process links this approach 
to some of the open systems work already mentioned. Other papers to recently 
use an entropy approach include the work by Lahiri [85]. 


h 


S [a^class] 


^finah^final 
^init »^init 
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The links between classical chaos and information theory are well known. The 
ideas can be applied to quantum mechanics too. Some recent work in this field 
has been done by Inoue et. al. [68,69]. See the references therein for historical 
details. 

Stochastic wehs (kicked systems) 

The pioneering work of Zaslavskii (see [24,116] for results and background) in¬ 
vestigating the effect of classical phase space structures on quantum dynamics for 
kicked systems has been fruitful. The ideas of the quantum suppression of chaos 
are clearly seen in many of these works through an analysis of the diffusive be¬ 
haviour. See the work by Berman et. al, Borgonovi and Rebuzzini, Chemikov et. 
al., Daly and Hefernan, Dana, Frasca, Korsch et. al, Sikri and Narchal, Torres- 
Vegaet. a/., and Zaslavskii et. al. [7,14,24,30-32,45,82,103,104,111,117,118] 
for just a few of the results obtained using these ideas. 

2.3.5 Spectral analysis of operators 

Finally, as alluded to throughout this review, the analysis of the spectrum of certain 
operators can be related to the dynamical properties of a quantum mechanical 
system. I mentioned this when discussing quantum billiards and random matrix 
theory but it also plays a central role in the examination of kicked systems. 

The spectral analysis of operators is a rich mathematical field in its own right 
and makes numerous claims relevant to the dynamics of quantum systems. The 
mathematical link between the spectral properties of operators and quantum dy¬ 
namics is the motivation for the work in this thesis and, accordingly. Chapter 5 is 
devoted to a fuller exploration of this link. 
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2.4 Quantum chaos—some established results 

Having reviewed some of the broad researeh areas in quantum ehaos, I now briefly 
present a number of established results from analytie examinations, numerieal 
simulation and experiment that have so far been established. They eut aeross view 
points in the field of quantum chaos. 

A common property of the time evolution of quantum systems is that for a 
short time the classical and quantum evolutions correspond. This correspondence 
is measured by, for example, the energy as a function of time. This is a generic 
property, and can be attributed to the time it takes for the quantum system to “be¬ 
come aware” of the finite dimensionality of its phase space. After such a time, 
the correspondence is lost and, while the classical system’s energy continues to 
increase (either in a diffusive way, corresponding to chaotic motion, or in a ballis¬ 
tic fashion, corresponding to resonant energy growth), the quantum system shows 
recurrences. A great amount of work has been done both analytically and numer¬ 
ically in identifying these timescales. See for example [7,58,70,76]. 

Related to this is the heuristic link between quantum recurrence and the prob¬ 
lem of conduction of electrons in a random lattice—the phenomena of Anderson 
localisation. Again, after a certain time, the classical and quantum evolutions 
diverge. See [22,43,44,51] for an explanation of the relation between quantum 
recurrence and Anderson localisation and also the contrast between classical and 
quantum time evolutions. 

The quantum recurrence results are so pervasive that they have lead to a new 
concept, already mentioned, the quantum suppression of chaos [12]. The concept 
encompasses all these results and reflects the fact that the quantum equivalent of 
many classically chaotic systems seem to be more well behaved and thus non- 
chaotic. In a general sense, this behaviour is attributed to the interference effects 
in the quantum evolution conspiring to suppress dynamical spreading in the wave 
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packets. This take on the results is beautifully explained by Haake [58]. He also 
shows how the introduction of tiny environmental interaetions destroys these in¬ 
terferences effeets, leading to rather chaotic-like behaviour for the open quantum 
systems. 


2.5 Summary 

I have diseussed eharaeterisations of ehaos in both elassical and quantum the¬ 
ory. While many of the simple ideas from classieal meehanies cannot be direetly 
translated into the language of quantum mechanics, there are ways around these 
problems, e.g., both the de Broglie Bohm theory and the Loschmidt echo approach 
allow concepts from classieal ehaos to be brought to the quantum theory. 

The “kieked” systems were seen to be of fundamental interest in investigations 
into quantum ehaos, as was the broad field of speetral analysis. Chapter 3 now 
introduees the basie tools required to eonsider the quantum evolution of kicked 
systems. 
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The Floquet Operator 
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The kicked systems are an excellent “testing ground” for chaos in both clas¬ 
sical and quantum mechanics. In the quantum case, the time evolution of the 
system, as governed by the Schrodinger equation, has a particularly elegant form, 
allowing a stroboscopic analysis of the system to be made. The mathematical tool 
for this is called the Floquet operator and is simply the time-evolution operator 
over a single kick period. It also goes by the name of quasi-energy operator or 
monodromy operator in the literature. The central role that the Floquet operator 
plays in the work presented in this thesis warrants a detailed introduction, the topic 
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of this chapter. For more information and some simple examples of the Floquet 
operator, see [58]. 


3.1 Time evolution in quantum mechanics 


The dynamical evolution of a non-relativistie elosed quantum system |'0(f)) is 
governed by the famous Schrodinger equation, 

i d 


hdt 




(3.1) 


With this notation, we have a rather simple looking partial differential equation. 
Of eourse, this is only an illusion. The vast complexities of quantum meehanies 
remain hidden. 

The formal solution to (3.1) is given by 

W)) 


exp { J dt' Hit'] 


IV'(f)) = 


where the “+” subseript indicates that we must ensure that the time ordering is 
done correctly. With correct time ordering, it becomes evident that the relation 


U{t,0) = U{t, s)U{s,0) 


holds and, due to the unitary nature of the time evolution. 


(3.2) 


U{t,s) = U{t,0)W{s,0). 


Apart from for the simplest of systems such as the harmonic oscillator, square 
well or the hydrogen atom (and minor variants on them), the Sehrodinger equation 
remains immune to analytic solution. Numerical studies get us some way further, 
but even then, all but the most simple of experimental situations remain intractable 
to detailed analysis. 
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The Hamiltonian H{t) can, in general, be explicitly time dependant, as in¬ 
dicated here. Solutions are just that bit more difficult. In some idealised cases 
however, there is an effective way to make progress and it turns out to be very 
useful in the study of potentially chaotic systems. As discussed in Section 2.1, 
the study of classical systems which experience a periodic (in time) perturbation 
is of great interest as the motion is often chaotic. If the Hamiltonian is given by 
some easily solvable system, say a harmonic oscillator or a quantum spinning top 
or rotor which is then periodically perturbed, the exact solution to the Schrodinger 
equation can be written down. If 

OO 

H{t) = Ho + (3.3) 

n=0 

where Hq is the time-independent Hamiltonian of the simple system, W is an op¬ 
erator describing the perturbation and A is a strength parameter, then the solution 
to the Schrodinger equation at times nT is given by 

l^(nT)) = V^\m), 

where I have introduced the so called Floquet operator V. V h simply the time- 
evolution operator for one period T, 


V = U{T). 


Using (3.2), the time evolution over one period is 


V = exp 


1 


r-T-e 


H{t) dt 


exp 


1 


fT-re 




H{t) dt 


The first factor, for time e < t < T — e, is trivial. The delta function in (3.3) is 
not acting, H{t) = Hq is independent of time and the system evolves freely via 
the time-evolution operator 


Uo = 
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The second factor is over an infinitesimal time period 2e when the delta function 
kick is acting. Over this infinitesimally short period of time the influence from Hq 
is zero. The system instantaneously evolves via the operator 

^-iXWjh 


Recombining these two parts of the evolution, the operator describing the evo¬ 
lution of the system from just before one kick to just before the next is seen to 
be 


Y ^ ^-iHoTin^-awin 


At times in this work, I will have need to consider the time-evolution operator 
taking us from just after one kick to just after the next. By a similar argument, this 
is given by 

y! ^ ^-iXWIh^-iHoTIh^ 


At any rate, they are essentially equivalent—in a given context, one may be more 
convenient than another. When used, it will always be made clear how the Floquet 
operator has been defined. 

Once the Floquet operator has been obtained, numerical investigations of the 
system become far more tractable. A stroboscopic picture of the evolution of com¬ 
plex systems can be obtained and investigated. It is also hoped that the idealisation 
of the delta function pulse is useful when it comes to predicting the behaviour of 
experimental situations where, for example, a system may be perturbed by a peri¬ 
odic stream of short powerful laser pulses. While each laser pulse clearly interacts 
with the system over some finite time, the effect on the dynamics should be mod¬ 
elled well by the Floquet operator type “kicks” discussed here. 

The Floquet operator also proves useful in analytic work. Specifically, the 
spectrum of the Floquet operator is of great use in discussing the dynamics of 
a given quantum mechanical system. Characterisations of the spectrum for par- 
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ticular classes of Hamiltonian systems may be of great value when it eomes to 
predieting how systems behave. 


3.2 An example: the kicked top 


Here I introduce the simple example of the quantum kieked top to demonstrate the 
usefulness of the Floquet operator. The analysis is covered in signifieantly greater 
detail in the book by Haake [58]. 

Both the unperturbed Hamiltonian, Hq, and perturbation operator, W, from 
(3.3) are polynomial funetions of the total spin of the top, J. Conservation of 
j +1)> where I have set ^ = 1 for eonvenienee, implies that the kieked top 
exists in a finite-dimensional Hilbert space. Choosing, somewhat arbitrarily (but 
for eonvenienee, see [58]), Hq oc and W (x. the Floquet operator is simply 
(ci and C 2 are the proportionality constants) 

Y — Q-iciJxT ^-i\c2Jl/2j 


where the l/2j factor in the instantaneous evolution term is required for the j —> 
oo elassieal limit to make sense. Again, see [58] for a diseussion of this. 

To proceed, we must speeialise to a partieular spin system. The spin-1/2 sys¬ 
tem is trivial, so here I briefly deseribe the spin-1 system. For spin-1 systems, 
states in the Hilbert spaee are represented by 3-component spinors (eolumn vec¬ 
tors), and operators by 3 x 3 matriees whieh ean be deeomposed in terms of the 
Gell-Mann matriees Ao,...,8- We have {Ji)im = —ifieum [88], so 


J. 


^0 -i 0^ 
i 0 0 

^0 0 Oy 


— A2 
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and 


Jx 


^0 0 0 ^ 
0 0 -i 

^^0 i 0 y 

Ay- 


Noting that 

= 2/3/ + Va/SAg 

the Floquet operator is 

'Y _ g-iciArTg—iAc2(2/3/+\/3/3A8)/2 
— g—iciArTg—ic3g—ic4A8 


where the new eonstants cg 4 have absorbed the various numerieal faetors and the 
kieking strength A. In a numerieal study using this operator, the strength is varied 
by adjusting C 4 appropriately. Note that the term is simply a global phase 
and henee has no effeet on the dynamies. 

The final step to be taken is to rewrite the Floquet operator direetly as a 3 x 3 
unitary operator. This ean be done either analytieally or numerieally. 

Now that we have an appropriate expression for the Floquet operator, the dy¬ 
namies of a quantum kieked top ean now be examined. An initial state, |^o) is 
ehosen 

(A 


1 ^ 0 ) 


^2 


vy 

and evolved using the matrix representation of the Floquet operator. The resulting 
states, lAinT)) are given by 


mnT)) = V^lAo) 
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and can be trivially generated numerically for further analysis. 

As 1/ is a 3 X 3 unitary matrix in this case, its spectrum is simply a point 
spectrum of three eigenvalues. 

3.3 Summary 

The unitary Floquet operator just introduced provides one with a stroboscopic 
view of the time evolution of periodically kicked systems. This was demonstrated 
through the simple example of the spin-1 quantum kicked top. The Floquet oper¬ 
ator is the key tool in the work that follows, but in order to understand the work 
a number of mathematical concepts must first be introduced. These concepts are 
the topic of Chapter 4. 
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As mentioned in Chapter 3, and to be detailed in Chapter 5, knowledge of the 
speetrum of the Floquet operator turns out to be important in understanding the 
dynamies of periodically kicked Hamiltonian systems. To make analytic progress 
in this field, we will need mathematical tools and concepts that are beyond those 
usually employed by physicists. In this chapter, I will provide a type of “tutorial” 
or, if you like, an overview of the fields of measure theory, functional and spectral 
analysis, picking out those concepts that will be necessary when we come to the 
work contained in Chapter 6. I begin with a review of integration and measure 
theory and then proceed to discuss how the concepts of singular, absolute and 
point measures (to be defined) apply in the context of operators on the physicists’ 
Hilbert space. 

4.1 Measure theory 

The typical high school introduction to “integration” begins with the idea that 
the integral of a function gives us the “area under the graph”. A suitably well 
behaved function is split up into smaller and smaller intervals, and each interval 
approximated by a rectangle, whose height is given by the function evaluated 
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at the mid-point of the interval. The sum of the areas of the reetangles gives an 
estimation of the area under the funetion. By taking the limiting ease of an infinite 
number of infinitesimal intervals, the integral of the function is obtained. 

This idea of integration is abstracted, formalised and extended to become the 
Lebesgue integral of modem mathematics [94, 98, 99,101]. To go beyond this 
point, the concept of abstract measure theory is introduced. 

When considering the integral of a function f{x), one usually thinks simply 
of 

Jf{x) dx 

where dx is the measure or “size” of an interval. Each interval along the x-axis is 
given the same weight. An abstraction of this technique is to allow the weight or 
“size” given to each interval to be determined by a measure function a{x). Any 
positive, non-decreasing function, a{x) will do. a{x) need not even be continu¬ 
ous. The measure 

/ia([a, b]) = o:{b) — a{a) 

is formed and one obtains the Lebesgue-Steiltjes integral 

pb pb 

/ / dua = f da. 

J a J a 

If a{x) = X, we trivially obtain 

b]) = a{b) — a{a) = b — a 


and recover the usual Lebesgue integral 




fdx. 


In an intuitive sense, when integrating with respect to a measure function a{x), 
the contribution to the integral from an interval is proportional to the derivative 
da{x)/dx. For a given function a{x), if the slope of a{x) is zero for a particular 
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interval, then that interval will eontribute nothing, irrespeetive of the funetional 
value f{x) on that interval. 

4.1.1 Point measures 

If the measure funetion is 

a{x) = 9{x — x') 

x' 

where 9{x — x') is the usual step-function defined as 

0 if x < x', 
x> x' 

then, in rather loose, but intuitive notation 

da{x) = 5{x — x')dx. 

x' 

As we move along the x-axis, we obtain a contribution to the integral j f da only 
at the points x'. Each contributing “interval” has a Lebesgue measure of zero— 
reflected in the fact that in “ordinary integration” the contribution to an integral 
from a single point is zero, and thus, the occasional infinitely thin but high spike 
in a function does not matter. 

Our integral with the given measure above is now 

J f da = ^/(x) 

x' 

and we have converted the integral into a sum. 

Measures of this type are called “point measures” for obvious reasons. Com¬ 
pared to integrating a function with respect to the Lebesgue measure, a point mea¬ 
sure gives contributions to the integral from single discrete points, exactly where 
the Lebesgue measure does not contribute. We say that the Lebesgue measure and 
point measure are mutually singular. 
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4.1.2 Absolutely continuous measures 

The absolutely eontinuous measure is perhaps the simplest of measures apart from 
the standard Lebesgue measure. If a{x) is a smooth, eontinuous, everywhere 
differentiable funetion of x then we can write 


A 

dx 


a{x) = g{x) 


with g{x) continuous and well-behaved. The integral J f dais then 



and we essentially recover the standard Lebesgue integral, but now of the function 
f{x)g{x) rather than f{x). The absolutely continuous measure satisfies 


/io,([a, b]) = 0 if and only if 6 — a = 0. 


A simple example is to consider the function a{x) = (l/2)x^. Then 5 f(x) = x and 
each interval in the integral, instead of having a constant weight, is now weighted 
by its position on the x-axis. 

The absolutely continuous measures are often referred to simply as the con¬ 
tinuous measures but I will not do so here, for reasons that will become clear in 
the following sections. 

4.1.3 Singular continuous measures 

The point and absolutely continuous measures discussed above are both fairly 
straight-forward. They quantify concepts that we are already familiar with. How¬ 
ever, as is usually the case in mathematics, we can take things further. 

It turns out that one can construct measures, /!„, that contribute to integrals ex¬ 
actly where the Lebesgue measure contributes zero (rather like a point measure), 
but which are nevertheless continuous measures. This seems contradictory but an 
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example should clarify the idea. We start by defining a particular set, the Cantor 
set. 


Consider the subset S' of [0,1] given by 


S 



27’27 J 


U ... 


The Lebesgue measure of S is 



The complement of this set, C = [0,1] \ S, has Lebesgue measure zero and is 
known as the Cantor set. It contains an infinite number of points but has a size of 
zero—it is an uncountable set of (Lebesgue) measure 0. See Figure 4.1. 


0 1/3 2/3 1 


Figure 4.1: The Cantor set 


The Cantor function, a{x), is defined by setting 

a{x) = 1/2 on (1/3, 2/3), 
a{x) = 1/4 on (1/9, 2/9), 
a{x) =3/4 on (7/9, 8/9)... 

a{x) becomes a continuous function by continuing this idea all the way to a{x) 
on [0,1]. The Cantor function a{x) is a non-constant continuous function on [0,1] 
whose derivative exists almost everywhere (with respect to the Lebesgue measure) 
and is zero almost everywhere! 
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The Cantor function, aptly coined the “Devil’s staircase”, is shown in Fig¬ 
ure 4.2. The function has zero slope almost everywhere, but still manages to rise 
from 0 to 1 across the finite interval [ 0 , 1 ] without overjumping by a finite amount. 



Figure 4.2: The Cantor function 


The measure, /Tq, formed using the Cantor function is continuous. That is, 
= 0 for any set {p} with only a single point, pa contributes to the 
integral only on the set C, exactly the set on which the normal Lebesgue measure 
contributes nothing. When integrating with this measure, almost every interval 
in [ 0 , 1 ] contributes nothing to the integral as the slope of a{x) is zero almost 
everywhere. 

The Cantor measure is an example of a continuous measure, but one which is 
mutually singular to the Lebesgue measure. It satisfies 

/ia([a, 6 ]) = 0 if and only if 6 — a 7 ^ 0 . 


Such measures are known as singularly continuous measures. 
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4.1.4 Summary 

I have introduced three types of measure—the point, absolutely continuous and 
singularly continuous measures. It turns out that an arbitrary measure can always 
be decomposed into three parts, 

/^PP /^ac /^sc- (4.1) 

This decomposition is unique and, importantly, the three pieces are mutually sin¬ 
gular. It should be noted that there are other ways to decompose the spectrum into 
mutually singular parts. While important, in this work I will not have need for 
such characterisations. The interested reader is referred to [94]. 

4.2 The tools of functional analysis 

Having introduced the notion of a generalised measure of a set on the real line, I 
now provide an overview of the field of functional analysis, through the introduc¬ 
tion of the Hilbert space and operators on the Hilbert space. The abstract mea¬ 
sure theory just presented is brought into the picture when discussing the spectral 
properties of operators on the Hilbert space. It is via this route that the power of 
functional and spectral analysis enters into mathematical physics. Classifications 
of the spectrum of operators (the eigenvalues) corresponds to finding the relevant 
measure to describe the eigenvalues. 

Here I only provide a brief overview, defining the pertinent concepts for my 
work and introducing notation. For a mathematically rigorous, step by step deriva¬ 
tion of the concepts discussed. Chapters 1 and 2 in [94] provide an excellent intro¬ 
duction. Much of the chosen notation is inherited from the work of Howland [66]. 
When, as will sometimes be the case, I need to depart from this notation, I will so 
indicate. 
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4.2.1 Hilbert spaces 

A vector space which is complete, i.e., one in which all Cauchy sequences con¬ 
verge to an element of the space, is called a Hilbert space. The elements of the 
space are the vectors which physicists use to represent states of a quantum me¬ 
chanical system. Throughout this work, I will always denote Hilbert spaces by H 
or /C. Subspaces of a given Hilbert space will often be referred to as S. A Hilbert 
space is said to be separable if and only if it has a countable orthonormal basis. 
I will only ever consider separable Hilbert spaces in this work. Elements of the 
space will generally be referenced by either |'0), 10), or in a mathematical, rather 
than physical context, may often be denoted by x or y. 

A subset S of H may or may not be complete. The closure of S, S is obtained 
by adding to S all the limit points of sequences of elements of S. By closing S, a 
complete subspace is obtained. A set S is said to be dense mHif S = H. 

The inner product of two vectors x,y E His {x, y), and the norm of a vector 
X is ||x|| = {x, 

4.2.2 Operators 

An operator A : H fC acts on elements of the Hilbert space H and returns 
elements in the Hilbert space /C. It is often the case that K, is either a subspace 
of, or in fact is, H. Operators on a Hilbert space will always be referenced by 
uppercase Latin characters. 

For an operator A : TY /C we define 

• the domain D{A) ■, the vectors x E H for which Ax is defined, 

• the range R{A) = {y E IC : y = Ax for some x E H}, 

• the kernel ker(A) = {x E H : Ax = 0}, and 
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• the operator norm IIAII = sup^g^(^), {Pa;||}. 

An operator is said to be densely defined if R{A) is dense in H. 

In this work I will have need to eonsider the family of all operators that ean aet 
on a Hilbert spaee H and produee an element of K,. This space, denoted /C) 
turns out to be a Banach space (a complete normed linear space). The Banach 
space is a generalisation of a Hilbert space, and thus some of the properties of 
Hilbert spaces discussed here do not apply to the space of operators. However, 
many do. Banach spaces will often be denoted by X and Y but context will avoid 
confusion with operators. See Chapter 3 of [94] for an introduction to Banach 
spaces. 

4.2.3 Invariance and reducing operators 

A subspace, S, of R. is called invariant for an operator A if, for all x E S, Ax E S. 
That is, action with A on 5 does not take us out of S. Further to this definition, 
a set 5 G 7^ is said to reduce an operator A if both S and its ortho-complement 
HQ S are invariant subspaces for A. 

4.2.4 Cyclicity 

A vector is cyclic for an operator A if and only if finite linear combinations of 
elements of are dense in H. This motivates the definition that a set S is 

cyclic for H if and only if the smallest closed reducing subspace of H containing 
S is H. This deserves some explanation. The existence of a cyclic vector means 
that by acting with the operator A and taking linear combinations of the results, the 
whole Hilbert space H can be explored. Now consider the reducing subspace S. 
Action on elements of S with A always leaves us within S. If repeated operations 
of this fashion end up reaching all vectors in the full Hilbert space H (i.e., there 
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are cyclic vectors in S) then the smallest reducing subspace of H containing S is 
H itself. Hence, the definition above for a set S to be cyclic corresponds to the 
existence of at least one cyclic vector in S. 

4.2.5 Limits 

If An is a sequence of operators, s-lim An (also An A) denotes the strong limit, 
defined by || [An — A)g\\^ 0 for all g w-lim An (also An A) denotes the 
weak limit, defined by | {Ang, f) — {Ag, f) \ — »■ 0 for all gf, / G H. By the Schwartz 
inequality, the weak limit exists if the condition is satisfied for f = g. I will also 
have need for the norm limit of an operator. An ^ A, defined by ||y4„ — y4||— 0. 
As the names suggest, 

An^ A^An^A^ An ^ A 

but the converse need not be true. A trivial application of the Schwartz inequality 
demonstrates this. 


4.3 The spectral theorem 

I refer the reader to Chapters 4, 6 and 7 in [94] for a full discussion of the following 
results. Here I will skip over much of the mathematical detail, only highlighting 
the pertinent definitions and theorems. 

The topic of spectral analysis is essentially an extension of the familiar linear 
algebra results on matrices of complex numbers to the action of operators on a 
Hilbert space—the spectrum is constructed in the same way that eigenvalues of a 
matrix are constructed. Many familiar mathematical expressions and definitions 
carry over, but great care must be taken to not incorrectly infer results based on 
those from finite dimensional linear algebra. The change from elements of C” 
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to arbitrary vectors in a (possibly infinite dimensional) Banach (or Hilbert) space 
allows for much more complex behaviour. 

I will start by defining the spectrum for an operator A : X X, mapping ele¬ 
ments of the Banach space X into the same space X. The space of such operators 
£(X, X) is itself a Banach space. Note that I generalise to consider the operator A 
acting on a Banach space, rather than a Hilbert space. It is normal practice in the 
literature to write £(X) for £(X, X) when the target space for A is also the initial 
space. With the spectrum defined, I will then jump ahead to the spectral theorem 
and the decomposition of the Hilbert space into point, absolutely continuous and 
singularly continuous reducing subspaces for a particular operator. 

4.3.1 The spectrum 

Consider A G £{X). For a G C, the operator 

{al-A)-^ 


is called the resolvent of A at a. 

The resolvent set, p{A), is defined by 

p{A) = {a G C : {al — A) is a bijection with a bounded inverse}. 

If a ^ p{A) then a is in the spectrum, a{A), of A. Essentially, the spectrum of 
A is the set of as for which the resolvent is not invertible. We now consider an 
element x of the Banach space X. If there exists an a G C such that 


Ax = ax 


then X is an eigenvector of A and a is the corresponding eigenvalue. For now, 
note that it is possible to have a G cr{A) without a being an eigenvector. See 
(p. 188, [94]). 
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4.3.2 Spectral measures 

I now specialise to operations on Hilbert spaces. As I will be concerned in this 
work with unitary operators, I state the following theorems in terms of unitary 
operators. The literature, however, generally introduces these concepts in terms 
of self-adjoint operators. As such, the references in this section refer to results on 
self-adjoint operators, but the results all pass through to unitary operators without 
much change. 

Fix A to be some unitary operator (typically it will be a time-evolution opera¬ 
tor) and choose ip ETi. As discussed in Section VII.2 in [94], one can show that 
there exists a unique measure, on the compact set a{A) with 

{'4’,f{A)ip)=[ (4.2) 

Ja{A) 

The left-hand side of the above equation is a typical inner product that a physicist 
needs to calculate. On the right-hand side is an integral in the complex plane, 
which we (hopefully!) have the tools to evaluate. The measure ij,^ is called the 
spectral measure for the vector ip. We now begin to see a link between physics 
and the mathematics presented earlier in Section 4.1 . The discussions have lead to 
a rigorous (although I have not really shown the details) derivation of the common 
physics practice of evaluating inner products via the insertion of a “complete set 
of states”. The flip-side is that beyond the standard discrete (point) and continuous 
states found in most physics books, we can now deal with arbitrary measures— 
point, absolutely continuous and singularly continuous measures. 

4.3.3 The spectral theorem 

First, I briefly return to the notion of cyclicity introduced earlier. Given an opera¬ 
tor A on a separable Hilbert space H, one can always find a decomposition of H 
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into reducing subspaces 

N 

^ = 0Hn 

n=l 

such that for each subspace Hn there exists a G Hn that is cyclic for A re¬ 
stricted to operating on Hn- Note that in the above, N can be finite or countably 
infinite. See (Lemma 2, p. 226, [94]). 

Defining T to be the unit circle in the complex plane, the spectral theorem 
states that for a unitary operator. A, on H, there exist measures {fJ-n}n=i = 
1 , 2 ,... or oo) on a{A) and a unitary operator 

N 

U :H^^L\T,dfin) 

n=l 

such that 

{UAU~^'ijj)^ (a) = aV'n(a) 

where -0 is a iV component vector ('0i(a),..., ipN{ot)) G ^n=i d/j^n)- Es¬ 
sentially, there exists a unitary operator U that converts operation with A onH 
into a linear combination of multiplications by a on elements of the complex unit 
circle. 

After introducing a few more concepts, I will return to another form of the 
spectral theorem which I will actually use for the rest of this work. It turns out 
that it is also the most familiar to physicists and provides a strong conceptual link 
between the mathematics here and the usual treatment of quantum mechanics. 

As discussed in Section 4.1, /ipp^ac.sc are mutually singular, and thus 

L^(R, dfi) = L^(M, d/ipp) © L^(M, dfi^c) © L^(M, dfisc) 

where I have introduced the space (M), the completion of the space C (M) of real 

continuous functions.* The spectral theorem says that operation with an operator 
'The completion implies the existence of a metric (or “distance” function). It is simply \\x-y\\ 
for x,y G H. It is normal practice to drop the measure, dji, when there can be no confusion, or 
when it is not necessary for an understanding of the concept being discussed. 
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A on 7 -^ is always equivalent to multiplication operations on the space L'^{a{A)) 
for some appropriate measure /r. Thus, the decomposition of L^(M, dfi) into point, 
absolutely continuous and singularly continuous parts allows one to conclude that 
for an operator A acting on H, one can always write 

H = Hpp © Tfac ® Tdsc 

where 

Hpp = point} = {an : is an eigenvalue of A}, 

Tfac = ii’lfJ'ip is absolutely continuous}, 

Tfsc = is singularly continuous}. 

Each subspace reduces the operator A. Introducing the notation A \ TYx for the 
restriction of A acting only on the subspace we may conclude that A \ Hpp 
has a complete set of eigenvectors while A \ Hac and A \ have only absolutely 
continuous and singularly continuous spectral measures respectively. 

The pure point, absolutely continuous and singularly continuous spectrum for 
an operator are now defined by 

(Tx(A) = a{A \ Tfx). 

It is possible to combine the sets in various ways to produce the singular (point 
plus singular continuous) spectrum and the continuous (absolutely continuous 
plus singularly continuous) spectrum 

(Ts(A) = crpp(A) U crsc(A), 

^cont(^) Crac(^) U (Tsc(^). 

As already alluded to, the form of the spectral theorem above is still not quite what 
is required in this work. By introducing the concept of spectral projections I will 
now rewrite the spectral theorem in a way that is very familiar to the physicist. 
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4.3.4 Spectral projections 

Section VII.3 in [94] covers the following rigorously as does Chapter VII in [98]. 

Spectral projections are, intuitively speaking, the building blocks for “operator 
measures”. Just as we considered the integral 



and characterised (decomposed) the measure function a{x) in terms of point and 
continuous measures, here I wish to decompose an operator A on the Hilbert space 
and write it as 



A 


where E(x) is an “operator measure”. 

The characteristic function Xs{x) for S' C M and x G M is defined by 



1 if a; G S', 


0 if a; ^ S'. 


The characteristic function provides the foundation for the step-function already 
introduced in Section 4.1 and is essential in a rigorous development of integration. 

The notion of the characteristic function is now extended to operators. To 
extend the definition, we must first understand how to form functions of operators. 

A function of an operator is, in general, defined through a limiting process 
of polynomials. Just as for functions on the real line, a general function of an 
operator can be considered as the limit of an appropriate sum of powers of the 
operator. The characteristic function is just one example. See [94, 98] for the 
details and the justification for why this formalism is self consistent, i.e., why it 
makes sense to talk of functions of an operator in this way. Note that this fact is 
essentially what allowed us to write down (4.2) earlier. 
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For a unitary operator A and a Borel (think “measurable”) set G T, the 
spectral projection, Eq -.H ^ His defined by 

En = Xn{A). 

The spectral projections Eq are orthogonal projections as E^ = Eq. A family of 
spectral projections, {Eq}, has a number of important and rather intuitive proper¬ 
ties that mirror the properties of ordinary measures. For a family of sets G T, 
the orthogonal projections have the following properties: 

• Each Eq is an orthogonal projection, 


• Eq — 0; Ej — I, 


• If with iln n ilm = 0 for all n ^ m, then 



Ecin Eci,n E* 




These properties are exactly those you would expect for a measure—a function 
that returns the “size” of a set. Any family of projections satisfying the above 
properties is called a projection-valued measure and is, rather confusingly, also 
denoted by Eq, where it is now understood that this is a family of spectral pro¬ 
jections, rather than a single spectral projection. As expected, it turns out that for 
any (j) e H, 

( 0 , En(p) 

is just an ordinary measure which one can integrate with respect to. The projec¬ 
tions, Eq, are self-adjoint^ and thus (0, E^j)) is real. 

^Actually, they are symmetric. An operator A is symmetric if, for /, 5 S , {Af, g) = (/, Ag). 
For bounded operators, as considered here, self-adjoint and symmetric are equivalent concepts and 
hence the use of “self-adjoint” is perfectly acceptable. 
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Armed with the speetral projeetions, we can now write an operator (introduc¬ 
ing the notation used for the rest of this work), A, simply as^ 



If A is unitary, then the eigenvalues are on the unit circle in the complex plane, 
and 



If A is self-adjoint, then 



A 


In the physicist’s language, the B(a) are the projections \4>n){4>n\ formed from 
the eigenvectors for the operator A. If the eigenvectors are discrete, A is 
expressed in terms of a sum. If the eigenvectors are continuous'^ then A is written 
in terms of a standard integral (with measure a the Lebesgue measure). 

4.4 Some results and discussion 

Having presented most of the concepts and definitions required for the rest of the 
work, I now summarise a few simple results that will be of use, and discuss some 
consequences of the theory so far presented. 

For any operator. A, the corresponding family F(a) form a general resolution 
of the identity 



^The subtle change of notation where I have replaced dE{a) with E{da) is a relic of the 


development of my work and the conflicting inherited notation from [ 94 ] and [66]. In the rest of 
this work, I generally use the notation E{da). 

'^We now know that it always eventuated that we had a discrete or absolutely continuous spec¬ 


trum in the standard physics examples. 
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Any function of an operator can always be written as 


f{A) = [ fia)E{da). 

Ja{A) 

As the E{da) are orthogonal projection operators, i.e., E‘^ = E, we have the very 
useful result that 


/ 


\f{a)fE{da) 


f{a)E{da) 


2 


The E{a), being operator measures, ean be deeomposed into their point, abso¬ 
lutely eontinuous and singularly eontinuous eomponents. Thus, we can decom¬ 
pose an operator A into three eomponents: App, A^c and Asc- We also form the 
singular and eontinuous parts of the operator A, 


A — A I A 

~ -^SC9 


A — 4-1-4 

-^cont ^ac ~ ^sc* 


Another minor ehange in notation due to the way in which this work developed is 
now introduced. As mentioned earlier in Seetion 4.2, in a mathematical context 
I usually use x and y to refer to elements of the Hilbert space H. The ordinary 
measure formed from a projection-valued measure E{a) for the operator A is now 
written as irix, rather than I define, for a vector x Edi, the measure 

mx{S) = {E{S)x,x) 

where S' is a Borel set in cr(A). Note that as E{S) is self-adjoint one is free to 
move it to the other side of the inner product. 

Having obtained a full charaeterisation of operators in terms of their deeom- 
position into point, absolutely continuous and singularly continuous parts, a final 
definition is now introduced. The operator A is pure point if and only if the eigen¬ 
vectors of A form a basis of Ti. That is, A is pure point if and only if TYcont = 0 
for the operator A. 
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4.5 Two examples—hydrogen and the harmonic os¬ 
cillator 

To clearly link the preeeding diseussion baek to familiar physios, oonsider the 
operator A to be the Hamiltonian, H, for hydrogen and oonsider the speotrum, 
aH,ofH. 

The bound states of hydrogen are a oountable number of isolated, diserete 
energies. Eaeh energy oorresponds to an eigenvalue of the system and the set of 
these points makes up the point energy speotrum. The positive energy seattering 
states form the oontinuous energy speotrum. Thus, the energy speotrum for the 
hydrogen system oonsists of two disjoint parts: the negative energy diserete (or 
“point”) speotrum, and the positive energy oontinuous speotrum. For hydrogen, 
ap{H) = {an] an ~ -13.6/n^ for n e N}, and acont{H) = (0, oo). 

The Hilbert spaoe splits into two subspaees Hpp and Ticont = Ti-ac- There is no 
singularly oontinuous eomponent to the Hamiltonian for hydrogen. This is typioal 
of most physies Hamiltonians. The singularly oontinuous eomponent of an oper¬ 
ator is rather abstraot and not oommonly oonsidered. In oertain oontexts however, 
we will see that it beeomes an essential tool in obtaining a better understanding of 
the dynamios of those systems. 

As another simple example, the harmonio oseillator quantum system has only 
diserete energy levels, and thus is said to be “pure point”. That is, the eigenveetors 
of the harmonio oseillator form a basis of the Hilbert Spaoe. 

The Hamiltonian of a typioal quantum meohanioal system does not posses 
a singularly oontinuous speotral eomponent and thus, singular oontinuity is not 
usually mentioned in texts on quantum meohanios. The physioal interpretation of 
a singularly oontinuous eomponent to the energy speotrum is murky to say the 
least. 
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4.6 Summary 


I am now ready to move on and diseuss why the spectrum of the Floquet operator 
is an interesting object to study in the context of quantum chaos. Results on the 
characterisation of the Floquet operator spectrum in terms of properties of the 
perturbation to the base Hamiltonian can then be presented. 
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can argue that knowledge of the spectrum of the Floquet operator provides great 
insight into the dynamics of the system of interest. 

The links between spectral decomposition and dynamics is still an active field 
of research. Model cases provide concrete examples of the links, but the interplay 
between spectra and dynamics is far from fully understood. Here, I outline the 
main results and provide references to the eurrent researeh efforts. 

I first eonsider the case of time-independent systems where the speetrum of 
the Hamiltonian is important. The RAGE theorem is the starting point for the 
diseussion. Time-dependent systems are then eonsidered—in partieular, time- 
periodic systems whence the Floquet operator’s spectrum is of importance when 
considering the dynamics. The work of Yajima and Kitada [115] on RAGE-like 
theorems is discussed, as well as the work of Milek and Seba [90]. The important 
eontributions of Antoniou and Suehaneeki [2,3] are also noted. A good overview 
of the field is provided by Combes [27] and also East [86]. The introduetion in the 
paper by Enss and Veselic [40] provides the most physieally intuitive discussion 
of the RAGE theorem and is essential reading. 

Finally, I also comment on a point of ambiguity seen in a number of papers. 
I detail the problem in Section 5.3 and suggest a possible reason for why the 
confusion has survived for nearly twenty years. 

5.1 Time-independent systems 

For time-independent systems, characterisations of the dynamics is linked to the 
spectrum of the Hamiltonian through, in part, the RAGE theorem [97]. A large 
literature considers these systems and the link between dynamics and the spectrum 
of H. See [86] and the references therein. 

A quantum system is characterised by its energy eigenstates, requiring knowl- 
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edge of the spectrum of H. The time evolution of these states is then examined. 
Thus, one considers terms like where '0 is spectrally decomposed with 

respect to the operator H. 

The RAGE theorem simply states that states in the point subspace, Hpp{H), 
survive as time evolves, while states in the continuous subspace, 7-^cont(-ff), decay. 
Put another way, for systems with a potential, states in the point subspace are 
bound states and essentially remain in a bounded region of space as time evolves. 
States in the continuous subspace spread in space as time evolves—they are the 
scattering states. 

The RAGE theorem is as follows [97]. Eet 77 be a self-adjoint operator and C 
a compact operator. Then for all ^ G HcontiH), 



The square on the norm can be dropped by technical considerations [97]. 

To better understand the RAGE theorem, consider the case where (7 is a pro¬ 
jection, namely [40], 


C = F{\x\ < R). 


F{\x\ < R) is the multiplication operator in x-space with the characteristic func¬ 
tion of the ball of radius R. The RAGE theorem then says that for %jj G FLcota{H) 
and a wide class of potentials. 



for any R < oo. That is, as time evolves, the wave-function ^ has negligible 
components within an arbitrarily sized ball in space. It spreads to arbitrarily large 
spatial distances from its initial location. 

Conversely, for ^ G Hp{H), for all e > 0, there exists an R{e) such that 


sup ||F(|x| > R) exp{—itF[)'ijj\\< e. 
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That is, there is some ball of radius R sueh that the wave-function has negligible 
components outside that ball for all time. The wave-function remains localised in 
space when ^|; G Hp{H). 

The RAGE theorem provides a course grained overview of the possible dy¬ 
namics. For 'll; G 7-fcont(-f^) the rate of decay is further dependant upon the char¬ 
acterisation of the continuous spectrum. Singularly continuous states are typi¬ 
cally expected to provide the “chaotic” behaviour, sitting somewhere between the 
ballistic (continuously accelerated) absolutely continuous states and the recurrent 
point states. 

5.2 Time-dependent systems 

First and foremost, I refer here to the work by Yajima and Kitada [115]. For 
systems with a time-periodic Hamiltonian (e.g., kicked systems) a RAGE-like 
theorem exists, but now the spectrum of interest is that of the Floquet operator, 
rather than the Hamiltonian. 

The numerous discussions on the link between dynamics and the spectrum of 
H for time-independent systems then apply. 

Yajima and Kitada show that for time-periodic systems the Floquet operator 
spectrum determines the dynamics in the same way that the Hamiltonian spectrum 
determines the dynamics for time-independent systems. 

For time-dependent systems as described in Chapter 3, the following two re¬ 
sults [115] apply. The decomposition of the Hilbert space, H, is with respect to 
the Floquet operator, V. Just as in the work by Enss and Veselic [40], Yajima and 
Kitada’s work [115] is presented in the case where C is the projection operator 

C = F{\x\ < R). 
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ForV' G Hcont(^), 


N 


lim -y 11^1^-^112=0. 

Af^ioo N ^ 


n=0 


For '0 G 7Yp(i/), for all 0 < e < 1, there exists an R{e) sueh that 


inf \\CV-i:\\> {1 - e) U\\ . 

n 

Simply put, for '0 G 7-fp(l/), the norm does not deeay. After some number of 
kicks, n, the change in norm is arbitrarily small and the state, initially localised 
in space, remains localised. The state -0 is a bound state of the system. If -0 G 
’^cont(i^) the norm does decay and the state -0 is a scattering state of the system. 

The work of Hogg and Huberman [60, 61] takes the point spectrum part of 
the result one step further. They show that the full wave-function “reassembles” 
itself infinitely often when the system has a discrete Floquet operator spectrum. 
The energy also shows recurrent behaviour and the system is seen to be quasi- 
periodic. 

The recurrence of the energy means that for a system to be chaotic, charac¬ 
terised by a diffusive growth in energy, it must have a continuous Floquet spectral 
component. The exact requirements are a topic of ongoing research. An intuitive 
discussion is provided by Milek and Seba [90] but it does have its flaws. I discuss 
this in the next section. 


5.2.1 Milek and Seba’s work 

The singular continuous spectrum of the Floquet operator can, and does, exist in 
physical systems of interest. With an understanding that Milek and Seba meant to 
refer to Yajima and Kitada’s RAGE-like theorem in [115] rather than the RAGE 
theorem itself, the argument presented in Section II of their paper [90] shows that 
if a system possesses a singularly continuous quasi-energy spectrum then its en¬ 
ergy growth over time may be characteristic of a classically chaotic system. Thus, 
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establishing the existenee or otherwise of singular eontinuous speetra of the Flo¬ 
quet operator is seen to be of eentral importanee to the question of whether or not 
a quantum meehanieal system is ehaotie. It must be noted that the arguments pre¬ 
sented by Milek and Seba are aeknowledged to be anything but rigorous—a point 
clearly established by Antoniou and Suchanecki [2, 3] who split the singularly 
continuous spectrum (and of course Hilbert space) into two parts, 


_ -^Decaying ^ -^Non-decaying 


The Hilbert space decomposition is reducing. Vectors in the decaying subspace 
do not survive as f > cxd while those in the non-decaying subspace do survive. 
The singular continuous part can act like either a point spectrum or an absolutely 
continuous spectrum. On reflection, this is to be expected. It is true, as argued 
by Milek and Seba, that parts of the singularly continuous spectrum will display 
chaotic behaviour. But there are also parts of the singularly continuous spectrum 
which do not. The game is wide open in terms of the details and subtleties, but 
there is no doubt that singularly continuous states can manifest as “chaotic like” 
behaviour. 

With the caveat that the described behaviour is not guaranteed, I now detail 
the argument of Milek and Seba. It has been inherited from the argument in [22]. 

Consider the probability Pk,i{n) of exciting the state of Hq after n cy¬ 
cles/kicks to the state, 

= injwT- 

To obtain the probability amplitude Pk An) I decompose the Floquet operator V 
into its point, singularly continuous and absolutely continuous parts. As V is 
unitary we can write, for uj an eigenvalue, 

V\u:) = 
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e^^^{k\E{dw)\l) 


We have 

Pk,i = 

' 't{V) 

Using (4.1) gives three parts, 

Pk,i = J 

= P'^^{k\u;){u;\l)+ f A^^hA^)du;+ f c(/i,e(M)(^)- 

^ '^<^ac(V') '^crsc{V^) 

The transition k —> I will only occur if a state, labelled by oj, links them. That is, 
for the sum part above, there must be at least one state |ci;) such that {k\uj){uj\l) A 
0. This demonstrates the importance of the spectral nature of V. If U’s spectrum 
is discrete (only the sum above remains) then the |a;) states are localised (like S- 
function spikes) and hence they can connect only a few states of Uq- In this case 
one would expect recurrent behaviour. As already discussed, this is the essence 
of Hogg and Huberman’s work [60,61]. 

If the spectrum is absolutely continuous, the Riemann-Lebesgue lemma (The¬ 
orem IX.7, p. 10, [95]) shows that 


lim Pk.iiji) = 0 V/c, 1. 

n^oo 

The system is continuously accelerated and we obtain resonant energy growth. 
This is not chaotic behaviour as the system’s future evolution will remain pre¬ 
dictable. 

If the spectrum of V has a singularly continuous component the possibilities 
for the dynamics are more varied. The RAGE theorem, as presented in Reed and 
Simon’s book (Theorem XI.115, p. 341, [97]), allows for a very slow diffusive 
growth in the energy in the presence of a singularly continuous spectral compo¬ 
nent (p. 343 and Problem 149 on p. 403, [97]).' As pointed out by Antoniou 

'Note that I am referencing the time-independent theory, but as shown by Yajima and Kitada 
the results flow through to the time-dependent theory which is of interest here. 
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and Suchanecki [2], this behaviour is not guaranteed—a point which Milek and 
Seba [90] seem to have missed or at least glossed over. Anyway, it is certainly 
possible that 


M 




There is a very slow increase in the energy of the occupied state over time. As 
the first equation diverges we see that the system spends an infinite amount of 
time in the “lower” states. That is, the system energy does not grow resonantly, 
but continues to explore the lower energy states as time progresses. The second 
equation shows, however, that on average the system will eventually escape any 
fixed chosen state—the system is not bounded in energy. 

This slow energy growth (diffusive) behaviour is typical of that seen in classi¬ 
cally chaotic systems as mentioned in Chapter 2. 

It should be noted that if there is no singularly continuous spectrum, then the 
diffusive growth described is certainly not possible. The existence of a singularly 
continuous spectrum is a necessary condition for chaotic behaviour. It is not a 
sufficient condition. 

5.2.2 The quasi-energy self-adjoint operator, K 

I have shown in the previous chapters that time-periodic systems are characterised 
by the Floquet operator—the unitary time-evolution operator over one kick period. 

There is an alternative way to access information on the spectral properties 
and dynamics of such systems. Developed in papers by Howland [64,65,67], the 
self-adjoint quasi-energy operator, or Floquet Hamiltonian, 


K = -id/dt + Hit) 
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turns out to provide a different way to aeeess similar information to what I am 
seeking from the unitary Floquet operator V. The speetrum of K is direetly related 
to the spectrum of V, as clearly shown in (p. 808, [16]). 

K also plays a central role in the work, already discussed, of Yajima and Ki- 
tada [115]. To obtain the time-periodic equivalent of the RAGE theorem, Yajima 
and Kitada first introduce K, apply the original RAGE theorem results, and then 
convert them into the result on the spectrum of the unitary Eloquet operator. 

A significant amount of the work on time-dependent systems utilises the Elo¬ 
quet Hamiltonian, K, because, in some sense, working with V proved difficult. 
The large body of knowledge on self-adjoint operators provides a mature basis 
for proving theorems about K. As discussed in [16], the spectrum of K is eas¬ 
ily related algebraically to that of V, so results on the spectrum for K and V are 
equivalent. 

The trade off is that, especially for a physicist, E is a far more intuitive opera¬ 
tor than K. The abstraction involved in working with K must be balanced against 
the loss of a simple physical picture. 

5.3 A point of clarification 

It must be noted that the concept of spectrum is associated with a particular opera¬ 
tor. Typically, physicists talk of the energy spectrum, associated with the Hamilto¬ 
nian. However, all operators (e.g., Hamiltonian, Eloquet etc.) have a spectrum. A 
failure to realise, or at the very least to explain, this has lead to a number of poten¬ 
tially misleading papers (see for example [90], but they inherited their argument 
from [22]) which used results on the spectrum of the Hamiltonian in a discussion 
of the spectrum of the Eloquet operator. While it is probable that the authors are 
aware of the jump they have made, referring to the work of [115], rather than the 
original RAGE theorem, would be of significant benefit. 
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Similarly, papers such as that by Guameri [53] do not make it clear which 
operator they are referring to in discussions of the spectrum. The issue is one 
of language. When physicists refer to “the spectrum” in quantum mechanics, 
it is generally assumed they mean “the spectrum of the Hamiltonian”. That is, 
the spectrum has come to mean the “energy spectrum” in the language of quan¬ 
tum mechanics. This, however, is not the mathematical definition. The spectrum 
is associated with a particular operator, as clearly discussed in Chapter 4. This 
misunderstanding has lead to a number of the more “physical” papers in the liter¬ 
ature incorrectly drawing conclusions from the rigorous mathematical literature. 
It does however turn out that most of the conclusions arrived at are valid. The 
time-independent system theorems (such as the RAGE theorem) all have equiva¬ 
lent theorems in the time-dependent theory (e.g., the work of Yajima and Kitada, 
Hogg and Huberman). That no obvious numerical inconsistencies have arisen 
means that the subtle flaws and potential misunderstandings in the physical liter¬ 
ature have gone largely unnoticed for close to twenty years. 

I believe that it is of paramount importance that authors exercise great cau¬ 
tion when discussing the spectrum of an operator. It should be made clear which 
operator is being investigated, especially when it is not the Hamiltonian. 

The work of Yajima and Kitada [115] is, in this field, the key link due to the 
fundamental importance of the RAGE theorem for time-independent systems. Un¬ 
fortunately, it is scantly referenced outside the mathematical literature. Increasing 
awareness of this important work would be greatly beneficial. 

5.4 Summary 

It is with the application of the RAGE-like theorem in mind [115], that I under¬ 
took the following work on the analysis of the quasi-energy spectrum of the class 
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of Hamiltonians as defined by (3.3). The aforementioned work by Milek and 
Seba [90], utilising the rank-1 work of Combeseure, has shown the manifestation 
of singularly continuous spectra in numerical simulations of rank-1 kicked rotor 
quantum systems. The work now presented in Chapter 6 and Chapter 7 extends 
these results and provides a rigorous mathematical basis to numerical calculations 
on the time evolution of higher rank kicked quantum systems. 




Chapter 6 


Spectral Analysis OF Rank-N Perturbed 

Floquet Operators 
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This chapter constitutes the first part of the research undertaken in my PhD. 
The aim is to characterise the spectrum of the Floquet operator for kicked systems 
as defined by (3.3). The method used parallels the investigation into the spectrum 
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of the Hamiltonian itself undertaken by Howland [66] and relies on the mathe- 
matieal baekground presented in Chapter 4. This work is an extension of a result 
of Combescure [28]. The motivation for investigating the speetrum of the Floquet 
operator has been discussed in detail in Chapter 5. 

6.1 Outline and summary of results 

I will derive conditions on the time-periodic perturbations to the base Hamiltonian 
for the spectrum of the Floquet operator to remain pure point. Equation (3.3) is 
replaced with a more “technical”, but equivalent form, 

OO 

H{t) = Ho + (6.1) 

n=0 

where A is bounded, W is self adjoint and Hq has pure point (discrete) spectrum. 
In terms of (6.1), the Floquet operator is 

y ^ yA*WA/h^-iHoT/h_ .g_2) 


If A is a rank-1 perturbation, 

A = \'ip){'ip\ 

W = XI 

then I reproduce the work of Combescure [28]. The vector \^jJ) is a linear combi¬ 
nation of the orthonormal basis states, \(j)n), of the unperturbed Hamiltonian Hq 

OO 

IV') = ^an|0n)- (6.3) 

n=0 

Combescure showed that if ^ G Ii{Hq), that is if 

OO 

n=0 


< OC 


(6.4) 
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then the quasi-energy spectrum remains pure point for almost every perturbation 
strength A. I will generalise this result to all finite rank perturbations 

N N 

k=l k=l 

N 

W = Y,^k\'ipk){'ipk\ (6.5) 

k=l 

where G M and each vector \^pk) is a linear combination of the Hq basis states, 

I0n) j 

CXD 

^^(®fc)n|0n)- (6.6) 

n=0 

The states \^jJk) are orthogonal, 

{'4’k\^i) = Ski. (6.7) 

The basic result is that if each \%jjk) is in Ii{Hq), the spectrum of V will remain 

pure point for almost every perturbation strength. 

The perturbation for which I prove that the quasi-energy spectrum remains 
pure point is, in fact, more general than the finite rank perturbation presented 
above. The finite rank result is, however, the motivation for undertaking this work. 

Howland [66] showed that the Hamiltonian (6.1) has a pure point spectrum if 
the V'fcS are in li(IIo). Here I follow a similar argument, showing that the continu¬ 
ous part of the spectrum of V is empty, allowing one to conclude that the spectrum 
of V must be pure point. 

Before proceeding, it should be mentioned that there are alternative routes to 
results similar to those I present. As mentioned in Section 5.2, associated with 
the unitary Floquet operator V is the self-adjoint Floquet Hamiltonian K [64]. 
Utilising K allows the self-adjoint work of Howland [66] to be used directly. This 
was done by Howland himself [67]. As my work is a unitary equivalent to the 
work of Howland [66] the results obtained correspond to those determined in [67]. 
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The relationship between my work and Howland’s work [66,67] is similar to the 
relationship between the self-adjoint rank-1 work of Simon and Wolff [105] and 
the unitary rank-1 work of Combescure [28]. 

The techniques developed in this chapter provide new, general theorems ap¬ 
plicable to unitary operators and show that it is possible to develop the theory of 
the spectrum of time-evolution operators directly, without need for the techniques 
of [64] briefly mentioned earlier in Section 5.2. 

In Section 6.21 will present the main theorems of the chapter, concerned with 
establishing when systems of the form given by (6.1) maintain a pure point quasi¬ 
energy spectrum. Parallelling Howland’s paper [66] on self-adjoint perturbations 
of pure point Hamiltonians, the key ideas are those of (7-finiteness and the abso¬ 
lute continuity of the multiplication operator V. To establish the second of these 
concepts for the unitary case (remember that we are concerned with the spectral 
properties of the unitary time-evolution operator and not with the spectral prop¬ 
erties of the self adjoint Hamiltonian), I will require a modified version of the 
Putnam-Kato theorem [96]. This, and associated theorems are the topic of Sec¬ 
tion 6.3. Section 6.4 uses the results of Section 6.2 and Section 6.3 to give the 
final results, which are then discussed in Section 6.5. 


6.2 Spectral properties of the Floquet operator 

Let (7 be a unitary operator on 77 and let /C be an auxiliary Hilbert space. Define 
the closed operator A : H ^ IC, with dense domain D{A). For our purposes, A 
bounded on 77 is adequate. I work with a modification (multiplication by e*®) of 
the resolvent of U, 

F(0;U) = 


( 6 . 8 ) 
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and define for 6* E [0, 27r) and e > 0 the function : K, K,, 

f/, A) = AF*{e+-U)F{e+- U)A\ (6.9) 

where 9± = 9 ± ie. Let J be a subset of [0,27r). 


Definition 6.1 (U-FINITE) The operator A is U-finite if and only if the oper¬ 
ator Gfi9] U, A) has a bounded extension to JC, and 


G{9-, U, A) = s-lim Gfi9-, U, A) (6.10) 

40 

exists for a. e. 9 E J. 


We define the function 




- E- 

27r ' ^ 


,m(t+2e) 






\n=0 


1 1 - e-2^ 

27r 1 — 2e“^ cos(t) + e~‘^^' 


( 6 . 11 ) 


The limit as e —> 0 of 6^(t) is a series representation of the (5-function. The proof 
is based on showing that 


lim J dt = 0 

where g{t) = f{t) — /(O) and f{t) is bounded in (—vr, tt). We split the integral 
into three parts, We must assume that f(t) is continuous at f = 0 

(otherwise f f{t)5{t) dt is not well defined) so that 


\/p, 3.^ > 0 s.t. Vf, |f| < .^ we have \f{t) — /(0)| < p. 


We have 




4 


g{t)6ft)dt= / g{t)6fit)dt+ / g{t)6fit)dt+ / g{f)5fit)dt. (6.12) 

J-TT JJi 
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Consider the third term in (6.12). For ^ < t < tt, cost < cos^, so 
1 + — 2e“^ cos t > 1 + — 2e“^ cos 


Therefore 


g{t)S^{t) dt 


r , ,1 1 

~ ^ ^ 27r 1 — 2e“^ cos.^ + e~‘^^ 


dt. 


Sinee g{t) is also bounded in (—tt, tt), we have 

|5'(^)| < K for t G (— TT, tt) 

for some iC G M. Thus, 


git)6eit) dt 


< 


K 


1 — e 


-2e 


< 


27r 1 — 2e“^ cos ^ 

Ktt 1 — 

27r 1 — 2e“'^ cos ^ ' 


dt 


Taking e —> 0, the denominator 

1 — 2e“'^cos^ + —> 2(1 — cos^) ^ 0 

as ^ > 0. The numerator 1 — e > 0. Therefore, 


g{t)5^{t) dt 


0 as e ^ 0. 


Similarly, the first term in (6.12), 

f g{t)5e{t) dt —0 as e —> 0 

J — 77 

tends to zero. Now eonsider the seeond term in (6.12). We have 


g{t)5e{t) dt 




< j \gi't)\5,{t)dt 


as 6^{t) is a positive funetion. As f{t) is eontinuous at f = 0, \g{t)\ < g for all 
\t\ < so 

j \g{t)\5e{t) dt < g j 5e{t) dt 
<g[ 6e{t) dt. 
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where we have again used the positivity of (5e(t) to extend the limits of the integral. 
This integral is equal to one by ((3.792.1), p. 435, [49]) so we conclude that 


g(t)Seit) dt 




< g. 


which can be made as small as desired and thus has the limit zero. We have shown 
that 


lim / [/(f) - /(O)] 5^(t) dt = 0. 

<^10 J-IT 


Therefore, 


lim/ f{t)S^{t) dt = \im f{0)6^{t)dt+ [f{t)-f{0)]6^{t)dt 

"10 J-n "10 J-n J-n 

/ TT 

/(0)(5e(f) dt 

= /( 0 ) 

and we have proved that lin\jo Se{t) as defined in (6.11) is an appropriate series 
representation of the ^-function. 

Given (6.11) and the spectral decomposition of U, we may write 

5, (1 - (7e*®) = J E{de') 

= [ 6,{9 - d')E{de') 


2tt 


2tt 


E{de') 




n=0 


n=0 




n=0 


n=0 


= —[F(»+;(7) + f(«+;(7)-l] 
1 


= ^(l-e-^‘)F*(«+;F)F(»+;f/). 


(6.13) 

The existence of a non-trivial U -finite operator will have important consequences 
for the spectrum of the Floquet operator V. We introduce the set 

JV([/, A, J) = {9 e J : s-limGe(6*; [/, A) does not exist} 

40 
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of measure zero, whieh enters the theorem. I will often refer to this set simply as 
N during proofs. 

Theorem 6.2 If A is U-finite on J and R{A*) is cyclic for U, then 

(a) U has no absolutely continuous spectrum in J, and 

(b) the singular spectrum ofU in J is supported by N{U, A, J). 

Proof (a) Following Howland, note that the absolutely eontinuous speetral 
measure, my^{J), is the e —> 0 limit of (^6^ (l — f/e*®) y, y) for 6* G J.lfyeH is 
in R{A*), allowing one to write y = A*x for some x G /C, then 

(1 - UR^) y,y) 

= \im(Sjl-UR^) A*x,A*x) 

= lim^^^^ {AF*{e+-, U)F{e+-, U)A*x,x) 

ej,0 TT 

= lim-^^-^ {Ge{d'i U, A)x, x) = 0 

ej,0 71 

for a.e. 6 E J. The set y of veetors y for whieh = 0 is a elosed reducing 

subspace of H, and by construction contains the cyclic set R{A*) as a subset. 
Because y is invariant, finite linear combinations of action with f/” leaves us in 
y. Due to the cyclicity, these same linear combinations allow us to reach any 
y G H. Thus, the set 3^ of vectors y with mf{J) = 0 must be the whole Hilbert 
space H. So there is no absolutely continuous spectrum of U in J. 

(b) A theorem of de la Vallee Pousin ((9.6), p. 127, [101]) states that the singu¬ 
lar part of the spectrum of a function is supported on the set where the derivative 
is infinite. In our case, this corresponds to finding where my{d9) —> oc. We 
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calculate 

lim{4(l-!/e")9.9) = j 5(0 - 0')(E(d0’)y,y) 

= J 5(9- e')my(de’) 

= my{d6). 

Thus, rriy = m^y + is supported on the set where 

lim ((5e (l — 17e*®) I/, I/) = cx). (6.14) 

From the proof to part (a), if y = A*x then the limit (6.14) is zero for 6 E 
J, 9 ^ N, so rriy in J must be supported by N. The set of vectors y with 
m*(J n N^) = my{J ~ iV) = 0 is closed, invariant and contains R{A*), so 
must be H by the argument above. Thus, the singular spectrum of U is supported 
on the set N. □ 

Now define a new operator, Q{z) : 1C ^ IC, 

Q{z) = A{l-Uz)-^A*. 

Note that 

Q (e*®±) = AF{e±, U)A*. (6.15) 

Q{z) is clearly well defined for \z\ ^ 1. Proposition 6.3 shows that the definition 
can be extended to [z] = 1. 

Proposition 6.3 Let A be bounded. I/OeJ, but 9 ^ iV(f/, A, J), then 

(a) the operator Q (e*®) = A (l — f/e*®) ^ A* is bounded on 1C, and 

(b) one has s-limQ = n 

ej,0 ^ ^ ^ ^ 
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Proof, (a) Without loss of generality, take 9 = ^ {z = 1). By Theorem 6.2, 
^ CTpiU), so (1 — exists as a densely defined operator. As A is a 

bounded operator, it suffices to show that (1 — f7e*°) ^ A* is bounded. We have 

II (1 - U(f^+Y^ A*x||2 = (F(0+; U)A*x, F(0+; U)A*x) 

= {AF*{f)+-U)F{f)+-U)A*x,x) 

= (a(0; U, A)x, x) < C\x\^ (as 9 ^ N) (6.16) 

for some real constant C.\fy = A*x, noting U = j e~^^E{d9), we also have 

(^37^) (id^) 

In light of (6.16), e may safely be taken to zero to obtain 

/ ( i_V^e ) (r^) Wd)y^y) < C ||x||2< oo. (6.17) 
From (6.17), we have 

/ (n^) (r^) w»)%v) 

= ([1 — U]~^y, [1 — f/]“^r/) < C ||x||^< cxo (6.18) 

so y e D[{1 — U)~^]. Thus, Q(l) = A(1 — U)~^A* is defined on all /C and 
bounded. 

(b) For y ^ D {{1 — U)~^), we show that the difference between Q 
and Q tends to zero as e —> 0. Again, due to the boundedness of A, we need 
only show that 

II (^(1 - f/e*°+)“^ - (1 - (7)-^) A*a:|| 
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tends to zero. Consider 

-1 


{l-Ue-^Y y-{I-U)-^y 


(l-e-*V") "E{de)y- I (l-e-*®) ^ E{de)y 

^ ^ ' E(dO)y 


1 — 


1 — 

T 2 


1 - e-*®e-^ 


:,-ie 


{E{de)y,y). 


1 — e 

To show that this has a limit of zero, write the numerieal faetor in (6.19) as 

2 


-fN2 


(6.19) 


1 

1 

2 

e (1 — e ^) 

1 — 

1 — 


(1 - e-*®e-") (1 - e-Y 


(1 - e-) 


(1 — e*®) (1 — e“*®) \ 1 — 26“*^ cos 0 + e 


-2e 


Equation (6.19) now equals 

(1 




{E{de)y,y) 


( 6 . 20 ) 


1 — 2e“^ cos 9 + / (1 — e“*®) (1 — e*®) 

The first faetor is bounded and tends to zero for 0 ^ 0. The second factor is the 
measure from (6.17). Clearly, away from the origin, the integral tends to zero. 
About the origin, some care must be taken to show that there is no contribution to 
the integral. 

Using (6.11), we have 


(1 




(1 




1 — 2e“'^ cos 9 + 

On substitution into (6.20), we obtain 


y^2nd.{B). 


(1-e-f 


2ir 




dQe 9y{9) 


■d9. 


1 — 2(1 —cos 6*) 1 — d9 2(1 —cos 6*) 

The function 0e(6*) = J S^(9') d9' is the step function in the e 0 limit. For 
non-zero e it is positive, monotonic, increasing and bounded by unity. As 9 ^ N 
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I have also written my{d9) = gy{6)d6 for some well behaved positive funetion 
gy{9). By integration by parts (see p. 32, [73] for existenee eonditions, whieh are 
satisfied) we obtain 


(1 


n-2e 


-27r 




9y{0) 


2(1 — cos6') 


0.W 


d gy{0) 
dO 2{1 — cos 9) 


d9 


The first term within the eurly braees is clearly some finite value. The second term 
is less than 


d gy{9) 


d9 = 


9y{0) 


d9 2{1 — cos 9) [2(1 — cos6')J 

from the properties of the 0^ function mentioned above. As with the first term, it 
is clearly some finite value. Noting that 

40 1 - 


part (b) follows. 


□ 


Theorem 6.4 Let A be bounded and U-finite on J, with R{A*) cyclic for U. 
Let W be bounded and self-adjoint on K,, and define the Floquet operator, 

y ^ yA*WA/hjj_ 


Assume that for \z\ 1, Q(z) is compact, and that Q converges to 

Q in operator norm as e ^ Ofor a.e. 9 in J. Define the set 

M{U, A, J) = {9 G J : Q does not exist in norm}. 

Then 

(a) V has no absolutely continuous spectrum in J, and 

(b) the singular continuous part of the spectrum of V in J is supported by the 
setN{U,A, J)UM{U,A, J). 
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Proof, (a) For convenience, write the Floquet operator as 

V = (1 + A*ZA)[/, 

where Z is defined appropriately by requiring* eyi\){iA*WA/h) = 1 -\- A*ZA. 
Noting (6.8) and (6.15) allows one to define 

Qi (e'®) = AF{e- V)A* 

= A {1- A*. 

Consider some vector y' G H. Ay' = x G /C is defined for such y'. A*x = 
y" is some vector in TY. The cyclicity of R{A*) means that action with linear 
combinations of powers of U on y" allows one to obtain any y ^ H, the original 
y' being one of them. Thus, we have a construction of A~^, namely, operation with 
A* followed by the linear combination of powers of U. As y' was arbitrary, A~^ 
exists for all y E H. This allows one to introduce I = A~^A in what follows.^ 

We now proceed by use of the resolvent equation, 

Qi — Q 

= aI^ _^ 

\ 1 - 1 - 

= {ttW (11 + - IIW") 

= ^ , (A*ZAUey — 

(1 — 1 — j 

= Qi (ey ZAUA-^e'^Q {e^ . (6.21) 

^For the rank-N perturbation case where W = ^k\fk){fk \ and A = \ fk){'4’k\, 

we have Z = I]fci(exp(tAfe/?i) - l)\fk){fk\- 

^The particular choice of T as a projection in (6.5) does not have an inverse, but I will show in 

Section 6.4 that one can define a subspace of Ti on which R{A*) is cyclic, and apply this theorem. 
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Thus, briefly using L = ZAUA for clarity, we have 

LQi -LQ = LQ,LQ 
1 — LQ + LQi — LQiLQ = 1 

^ (l + LQi)(l-Lg) = 1 

^ 1 + e^^ZAUA-^Q^ (e*®) = 

[1 - e^^ZAUA-^Q{e^^)]~\ (6.22) 

Denote by N and M the sets N{U, A, J) and M((7, A, J). If 6^ G (J ~ A^) ~ M, 
i.e., 6* G J n n M^, and 1 — e^^ZAUA~^Q (e*®) is not invertible, then the 
compactness of —LQ (e*®) (which follows from the compactness of Q 
the norm convergence of Q and (Theorem VI. 12, [94])) allows one to 

use the Fredholm Alternative (Theorem VI.14, p. 201, [94]) to assert that 

3x G /C, s.t. [1 - e^^ZAUA-^Q (e*®)] x = 0. 


That is, there is some vector x ^ 1C which satisfies the equation 


X 


e^^ZAUA-^A (l - f/e*®) ^ A*x = 0. 


As 0 G J ~ A^, by Proposition 6.3 p = A*x e D (l — (7e*®) 

0 = (1 - A*X. 


(6.23) 

so define 0 as 

(6.24) 


0 is a well defined vector on H and we have 

X - e^^ZAUcj) = 0, 


which implies that 


X = e^^ZAUcj). 
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By (6.24), x 7 ^ 0 implies 0 7 ^ 0, so we have 

(1 - Ue^^) (j) = A*x = e^^A*ZAU(j), 
whenee (1 + A*ZA) (70 = e“*®0, 

or 1/0 = e-*V- (6.25) 

We eonelude that G (TpiV). 

The multiplieity of the eigenvalue is given by the dimension of the kernel of 
1 — ZAUA~^Q, whieh is finite by the eompaetness of Q and (Theorem 4.25, 
[99]). 

Therefore, if 6 >GJ~(A^UMU apiV)), which is a set of full Lebesgue 
measure,^ then the vector 

x{e) = [1 + e^^^+^^^ZAUA-^Qi x 

= [1 + L+Qi (e'®+)] X (6.26) 

must be bounded in norm as e —0 because we have just seen that if it is un¬ 
bounded we have an eigenvalue of the operator V. For y = A*x G R{A*), the 
absolutely continuous spectrum, m“, of V is the limit of 

(<5, (1 - y, y) = {A5, (l - A*x, x) . 

The aim is to show that this is zero for all y G 77. Define 

Fi(0) = (6.27) 

F{e) = {l-Uey~^ (6.28) 

and in a similar fashion to ( 6 . 21 ) and ( 6 . 22 ), we obtain 

Fi(0) = F{e) [1 + (1/ - U)F^F^{e)\ (6.29) 

^That the set M has measure zero is a consequence of Lemma 6.5 on page 79. 
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and 

(1 + (1/ - U)e^^Fi{e)) = (1 - (1/ - U)F^F{e))~^. 

Writing X = V — U, on substituting (6.29) into the expression for the (5-funetion 
(6.13) we obtain 

2716, (1 - = (1 - e-2^) F*{e^)Fr{e+) 

= [1 + F^+XFi{e+)]*27r6, (1 - UF^) [1 + F^+XFi{e+)] . 

Substitution of (6.15) and noting that 

X = V-U = {1 + A*ZA)U -U = A*ZAU 


gives 

(1 - VF^) ri* 

= ri [1 + XF^+Fi{e+)X 5, (1 - UFX [1 + XF^+Fi{e+)\ ri* 

= ri [1 + e-^^-Fl{e+)U*A*Z*A] 6, (l - UF^ [l + F^+A*ZAUFi{e+)\ ri* 
= [ri + e-^^-AFl{e+)A* (ri*)"^ U*A*Z*A] 6 , (l - UF^ 

X [A*+ F^+A*ZAUA-^AFi{9+)A*] 

= [l + e“*®-Qi(0+) (ri*)"^ U*A*Z*] A6, (l - UF^ 21* 

X [l + F^+ZAUA-^Q^{e+)] 

= [1 + LM0+)Y 215, (1 - UFX 21* [1 + L+Q,{e +)]. 

The absolutely eontinuous speetrum, of V is the e —> 0 limit of 


(215, (1 - 12e*®) A*x,x) 

= ([1 + L+Qi(0+)]* 215, (1 - UFX 21* [1 + T+Qi(0+)] x, x> 
= {A6,{1-UFX A*x{e),x{€)) 


e(l 


^{G,{e-,U,A)x{e),x{e)) 


TT 


(6.30) 
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which tends to zero as e —> 0 if both G^{9; U, A) and x{e) are bounded. 0^(9; U, A) 
is bounded as 9 e J ^ N and x{e) is bounded by (6.26). 

Part (a) follows since R{A*) cyclic for U implies that R{A*) is cyclic for V. 

(b) Let = N{V, A, J). We have just shown that 9 E J ^ {N UMUap{V)) 


implies that 




-{GF9-,V,A)x{e),x{e)) ^0 

71 

(6.31) 

and therefore 




{6,(1-ve") y,y)^0. 

(6.32) 


If we can infer the strong limit from this weak limit then we have established 
that 9 ^ A^i. We use the result that if Xn ^ x and ||x„||—>||x||, then Xn x 
( [6], p 244). Writing and G for G^{9-, V, A) and G{9-, V, A), and and F for 
F(6'+; V) and F{9] V), eonsider 

[lIGeXlP — ||Gx|P| 

= 1((G'^ - G‘^)x,x)\ 

= 1(21 {{F*F, - F*F) A*AF*F, + F*FA*A {F*F, - F*F)} A*x, x)l. 

If A, F^ and F are bounded operators, then if F*F^ — F*F tends to zero as e —> 0 
we ean conclude that the strong limit exists. A short calculation shows that 

F*F, - F*F = [(1 - e-2^) - (1 - e-^) F*F,F*F 

which trivially tends to zero as e —> 0 given the boundedness of F^ and F. Finally, 
A is bounded by assumption and (6.26) shows that Qi{9+) is a bounded operator 
as e —> 0 and thus both F^ and F are bounded. 

Moving on from (6.32), we have now established that C NUMUap{V) so 
Ni must have measure zero, again remembering that we need Lemma 6.5 below 
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to prove that M has measure zero. By Theorem 6.2, Ni supports the singular 
speetrum of V. That is, 

(iV[) = 0 

where the set iV[ is the eomplement of Ni. As the measure is positive and = 
m^c _j_ know that 

(iV[) = 0. 

Trivially, (iV U M) ~ crp{V) contains ~ crp(l^). Thus 
([iVi n (iV[ U ap{V)) 

= {N^)+m^‘^ MV)) 

= 0 + 0 
= 0 

as the (continuous) measure of single points is zero. 

The set A^UM r\ap{V)^ must support as Nir\ap{V)^ is a subset. Therefore 

m^^[N U M n apMf) = 0. 

This equals 

m"" {[N U Mf U ap{V)) = ([^ U Mf) + m"" MV)) 

so we conclude that the set iV U M supports the singular continuous part of the 
spectrum. □ 

Theorem 6.4 has shown us that V has an empty absolutely continuous compo¬ 
nent, and that the singular continuous component is supported by the set N U M, 
which is independent of A. We know that N has measure zero, and Lemma 6.5 
below shows us that M also has measure zero. This will allow us to apply Theo¬ 
rem 6.6 to show that the singular continuous spectrum of V is also empty. Thus, 
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with both the absolutely eontinuous and singularly eontinuous spectra empty, we 
can conclude that V must have pure point spectrum. 


Lemma 6.5 Let Q(z) be a trace class valued analytic function inside the com¬ 
plex unit circle, with \z\ < 1. Then for a. e. 9 

limQ = Q 

exists in Hilbert Schmidt norm. 


Proof We parallel the proof of de Branges theorem (see [34] and p. 149- 
150, [77]). Consider 


Q + Q* 

^[a 


+ \ AE{d9') 


= / yf* 


_ g-i(6»'-0)g—e 

2 (1 — cos(6'' — 6)) 


AE{de'). 


1 -I- — 2e“'^ cos(6*' — 6) 

The factor within the curly braces is greater than zero for all 6', 9 and thus we have 


Q + Q* (e*(®+-)) > 0 Ve > 0. 


Therefore, following de Branges, 

|det (1 + Q ( > det (1 + Q* (e‘<'’+“>) Q (e‘<'’+”>)) 

= n (1 + I“•>I7 


> 


1 . 


{an} are the eigenvalues of Q . From these two bounds we obtain 


Q (g40+*O) 


1 

det (1 -f Q (e*(®+*0)) 

H.S. 

det (1 -f Q (e*(®+*0)) 
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The definition of an analytie operator (p. 189, [94]) implies the analytieity of 
the eigenvalues, and thus the operations of taking the determinant and the Hilbert 
Sehmidt norm are analytic. Hence, both functions above are analytic and bounded 
within the complex unit circle (e > 0). Application of Patou’s theorem (p. 454, 
[36]) establishes the existence in the limit as e —> 0 and hence both functions exist 
on the boundary almost everywhere. Taking the quotient we establish the exis¬ 
tence of Q in the Hilbert Schmidt norm. □ 

Let (fl, /i) be a separable measure space, and 



a measurable family of unitary operators on H. We denote by 



the multiplication operator 


(Vm)(A) = 1/(A)m(A) 


on /i; H), where u{X) G L^(H, /r; H). 

A vector m(A) is an element of /r; H) if, for u{X) G H, 



It is important to note the difference between 1^(A) acting on H and V acting on 
/i; H). To obtain our goal of showing that for a.e. A, V (A) has a pure point 
spectrum, we must show that V is absolutely continuous as a function of A on the 
space /r; H). 

Theorem 6.6 is taken directly from [66]. The proof given is, apart from some 


small notational changes, identical to that in [66]. Due to a number of typograph¬ 
ical errors however, I have reproduced the proof here for reference and clarity. 
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Theorem 6.6 Let V be absolutely continuous on //; 7i), and assume that 
there is a fixed set S ofLebesgue measure zero which supports the singular contin¬ 
uous spectrum ofV (A) in the interval J for p-a.e. A. Then V (A) has no singular 
continuous spectrum in J for p-a.e. A. 


Proof For fixed x &TL, and any measurable subset F of Ll, let u{X) = Xr(A)x 
be a veetor in p; Ti). Then 

j \Eiy]x\^ p{dX) < j \Ex[S]xf p{dX) 

= J\E,[S]u{X)f p{dX) 

= J |E[5]M(A)|V(t^A) 

= ||E[%(A)|p = 0. 

/p \Ef^[J]xf p{dX) = 0 implies that \E{fJ]xf = 0 for /i-a.e. A. Thus 


ElfJ]x = 0 


for every x E H. □ 

The applieation of Theorem 6.6 relies on finding a fixed set S of measure zero 
whieh supports the singularly eontinuous speetrum. S = N U M is suffieient. 

I have now established all the basie requirements for V to be pure point, given 
U pure point. They are now eombined to produee the main theorem of the ehapter. 
There is still quite a lot of manipulation to satisfy the eondition V absolutely eon¬ 
tinuous on L^(M; Ti) of Theorem 6.6, and this will be the foeus for the remainder 
of Seetion 6.2 and Seetion 6.3. 



82 


Spectral analysis of rank-N perturbed Floquet operators 


Theorem 6.7 Let U and A satisfy the hypotheses of Theorem 6.4 and define for 
A G M 

y(A) = 

Then V (A) is pure point in J for a.e. A. 

Proof By Theorem 6.4, with W = XI, i^(A) has no absolutely eontinuous 
speetrum in J, and its singularly continuous spectrum is supported on the fixed 
set S = iVUM. Application of Lemma 6.5 shows that S is of measure zero. If we 
can show that V is absolutely continuous on L^(M; TL) then Theorem 6.6 applies 
and shows that the singular continuous spectrum is empty. I prove the absolute 
continuity of V in the following sections. 

As I have shown that both the absolutely continuous and singular continu¬ 
ous parts of the spectrum are empty, we conclude that 1^(A) is pure point for a.e. 
A G M. □ 

To show that V is absolutely continuous, I apply a modified version of the 
Putnam-Kato theorem which is proved in Section 6.3. The unitary Putnam-Kato 
theorem is: 

Theorem 6.11 Let V be unitary, and D a self-adjoint bounded operator. IfC = 
V\y* , D] > 0, then V is absolutely continuous on Hence, if is 

cyclic for V, then V is absolutely continuous on TL. 

I apply this theorem on the space L^(M; TL). A naive application to obtain the 
desired result is as follows. I slightly change notation and explicitly include the 
A dependence of W in the definition of V. If we choose D = —i{d/d\), with 
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V = then 



-iXA*WA 


-V*A*WA, 


so that for some u e L'^ 




[W*,3]u = (y*B-BW*)u 


-DV*m 


. d 


(V*m) =W*A*WAu. 
dX 


Therefore, 


C = V[V*,D] = A*WA. 


With VF = /, we obtain C = A*A > 0 and thus = R(A*) (see the proof 

to (Theorem VI.9, [94])) is eyelie for V. Hence, V is absolutely continuous and 
all the requirements of Theorem 6.7 are satisfied. 

The problem here is that D is not bounded, and boundedness of D is essential 
in the proof of the Putnam-Kato theorem. I use a similar technique as Howland 
[66] to overcome this issue. 

As the norm of A* A may be scaled arbitrarily, we can rewrite V, for real t, as 


V{t) = 


(6.33) 


for some real c > 0. 

Proposition 6.8 On L^(M;7-(), consider the unitary multiplication operator 
V, defined by 

Yu{t) = V{t)u{t) = 

and the bounded self-adjoint operator D = — arctan(p/2), where p = —id/dt. 
Then C = V[V*, D] is positive definite, and R (67^/^) is cyclic for N. Hence, the 
requirements of Theorem 6.7 are fully satisfied. 
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Proof. The operator D on L^(M; H) is convolution by the Fourier transform 
of — arctan(a:/2) [66], which is ((3), p. 87, [41]). This is a singular 

(principal value) integral operator, because arctan(p/2) does not vanish at infinity. 
Thus, for M(f) G 

j-oo ^-2\t-y\ 

DM(t) = inP / - uiy) dy 

7-00 t-y 


and 

[YfB]u{t)=iTTpJ g-2|t-y| ^ ^^'^ u{y)dy 

so 


Cu{t) =Y[Yf]D)]u{t) 

7-00 t - y 

Inserting expression (6.33) for V(t), we obtain 

roc -I _ pic(t—y)A*A 

Cu{t) = in / - u{y) dy 


(6.34) 


' —CXD 

poo 


= ITT 


t-y 

- 2 \t-v\ 1 - COS {A*Ac{t -y))-i sin {A*Ac{t - y)) 


t-y 


u{y) dy. 

(6.35) 


Note that this is no longer a singular integral. To show that C is positive, we must 
show that 

{u{t)Xu{t)) > 0 yu{t) e L‘^{R-,n). 

Note that the inner product on L^(M; 77) is given by 


{u{t),u\t)) 


u*{t)u'{t) dt. 


(6.36) 


The operator A is now written in terms of its spectral components. Note that here 
A decomposes A and bears no relation to the strength parameter used at other 
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stages in this chapter. When required for clarity, I write to identify the integral 
over the variable A, 

A = J XE{dX). 

A general vector u{t) may be written 

u{t) = J E{dX)u{t). 

Then 

f{A)u{t) = J fiX)E{dX)uit) 

which implies that we may rewrite (6.35) as 

roo P , 1 _ pic(t-j/)|A|2 

Cu{t) = m / dy - E{dX)u{y) 

J-oc J\ t-y 

= f dy f (j)x{t - y)E{dX)u{y) 

J — oc J A 

= j^E{dX)Cx{t) 


where 


C\{t)= / dy(j)x{t-y)u{y) 


and we have defined the new function 


Mt) = ivre-2|*lr^ (l - . 


By the convolution theorem, note that 


Cx{uj) = 0A(ct')fi(a;) 


where the indicates Fourier transform. 

Using this decomposition of u{t) and ParsevaTs theorem, we can now easily 
write down Cu{t)). I use (x, to indicate the inner product on the Hilbert 
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Space H, reserving (x, y) for the inner product on L^(M; H) as in (6.36). 

poo 

{u{t),Cu{t)) = I dt {u{t),Cu{t))^ 

dt (u{t), [ E{d\)C\{t) 


' —OO 

poo 


' —OO 

poo 


n 


’ —OO 

poo 


dt / E{dX) / dy - y)u{y) 


n 


n 


= J dt (^u{t), J^E{dX) j 
= j E{dX) j ^ (u{uj),^x{u:)u{u:)^^ 

= J^E{dX) J ^ \ux{uj)f 

It is clear that if 0 a (<^) is positive for all A then C will be positive. 

In the following calculation we will find the need to bound c|Ap. The restric¬ 
tion 0 < c| Ap < 1 will be employed. I argue that as ^4*^4 is a positive self-adjoint 
bounded operator we can restrict the integral over A to ([98], p. 262, 273) 


fM 


A* A = 


\X\^E{dX) = / \X\^E{dX) 


(6.37) 


' m—0 


where M is the least upper bound and m the greatest lower bound of A* A. The 
norm of A* A is given by max(|m|, |M|). Thus, by setting 

1 


c = 


IIAMII 

then each c| Ap is guaranteed to be less than unity. 

Proceeding, the Fourier transform, 0 a (<^), of 

(pxit) = [l — cos cf|Ap — i sin cf|Ap] 


(6.38) 


is now calculated. Split (6.38) into two parts, 

0Ai(f) = [l — coscf|A|^] , 

0A2(f) = vre“^'*0“^sincf|A|^. 


(6.39) 

(6.40) 
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The Fourier transform of (6.39) is 

/OO 

— coscf|Ap)e“*‘^* dt 

■OO 

poo 

= in / — cosct|Ap)e“*‘^* (if 

Jo 

poo 

Jo 

Using ([41], p. 157, (59)), and setting S = c|Ap/(2 + iuj), we obtain 

2 r \ 1 f ^ + 

■^"(‘") = Ti°g(TTs^, 

The logarithm of a eomplex number ean in general be written as 


log(z) = log(|z|) + i Arg 


so noting that | (1 + A^)/(I + S*"^) | = 1, we see that 

2 / \ ^ A / 1 + 5'^ 


= -vrArg (l + . 


With K = c|Ap, the real and imaginary parts of 1 + are 


3? (1 + S^) 
A (1 + 


{A + u^f + K^ (4-o;2) 
(4 + o;2)2 
—4fi;^a; 

(4 + o;2)2- 


With the restriction that 0 < k < 1, the real part is positive for all u and thus 

Arg(z) = arctan(Az/3?z). Thus, 


(j)xi{uj) = —Trarctan 


p(l + 5^) A 
1s(1 + S2)T 


arctan(z) is the principal part of Arctan( 2 ;), with range —n/2 < arctan( 2 ;) < 
n/2. The Fourier transform of (6.40) is similarly calculated, using ( [41], p. 152, 
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(16)), to be 


4>X2{^) = TT [arctanS' + arclans'*] 


= vr 


arc! an 


c|Ap 


+ arc!an 


2 + iu! J \2 — iuj 

Repeated application of the formula arc! an ( 2 : 1 )+arc! an ( 2 : 2 ) = arctan(^i+; 22 /l- 
Z 1 Z 2 ), valid when Z 1 Z 2 < 1 (true for 0 < k < 1), yields'* 


+ 0A2(ct') 

'n(o^,c|Ap] 


= vr arc! an 


d{(j, c|A|^ 


(6.41) 


where 


77,(0;, k ) = 4 k [(4 + + Ko; (4 + o;^) + (4 — o;^) — k^o;] 


(6.42) 


and 


d{u;, k) = (4 + - 2k^u;^ (4 + o;^) - 16 k^o; - (4 - . (6.43) 

We can easily confirm that for 0 < k < 1, n{uj, K)/d{uj, k) and hence is 

strictly positive by noting that there are four distinct regions of interest for oj, in 
which terms in n and d do not change sign. Table 6.1 shows these regions and the 
sign of each term in the region. Note that the global (positive and hence irrelevant) 
K factor from (6.42) is dropped from the numerator for the following discussion. 

■^This result is not valid for values of k larger than around 2, at which point the arctan addition 
formulas fail—this is a moot point however, as we may trivially restrict k as already explained. 
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n{uj, k) = 

{A + uj^f 

+ KUJ (4 + Uj"^)^ 

+^2 (4 - UJ^) 

— K^UJ 

u < -2 

+ve 

-ve 

-ve 

-tve 

-2 < U! <0 

+ve 

-ve 

-i-ve 

-tve 

0 < u <2 

+ve 

+ve 

-i-ve 

-ve 

UJ >2 

+ve 

+ve 

-ve 

-ve 

d{uj, k) = 

{A + uj'^f 

-2k^uj‘^ (4 + uj‘^) 

—IQk^uj 

(4 - uj^) 

UJ < -2 

+ve 

-ve 

-I-ve 

-tve 

-2 < a; < 0 

+ve 

-ve 

-tve 

-ve 

0 < a; < 2 

+ve 

-ve 

-ve 

-ve 

UJ >2 

+ve 

-ve 

-ve 

-tve 


Table 6.1: Sign of each term in the numerator n(ci;, k.) and the denominator d{uj, k) 
of (6.41). 


For each row in the table, we simply need to show that the terms add to produce 
a strictly positive number. First note that the first column for both the numerator 
and denominator is independent of k. To show the positivity of each row, set all 
positive K-dependent terms to zero and then take k = 1 for the negative terms 
to maximise their contribution. Expanding out terms, it is then trivially seen in 
all cases that the first column ((4 + for the numerator and (4 + for the 
denominator) dominates. Thus, no row is negative and we conclude that (px is 
positive definite. 

We have established that the Fourier transform of (px is positive definite for 
c|Ap < 1. As a visual aid. Figure 6.1 shows (pxi^)- The positivity for c|A| 2 < 1 
is clear. 
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Figure 6.1: Plot of (f)x{uj) = + (j)x 2 {to), the Fourier transform of 4>x{t) = 

(1 — cos(cf|Ap) — zsin(cf|Ap)). 0 a(<^) is strictly positive for all u 
when c|Ap < 1. 


Thus, C is strictly positive and V is absolutely continuous on As 

A* A is a factor of 1 — (i.e., A*A is a factor of C), = R{A*). Not¬ 

ing that R{A*) is cyclic for U and hence cyclic for V, we conclude that 
is cyclic for V. Thus, V is absolutely continuous on L^(M; H). □ 


I have now satisfied all the requirements of Theorem 6.7. 
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6.3 The unitary Putnam-Kato theorem 

In this section, I will prove a modified version of the Putnam-Kato theorem, as 
used in the preceding section. The theorems and proofs follow a similar argument 
to that of Reed and Simon (Theorem XIIL28, p. 157, [96]) and are motivated by 
the stroboscopic nature of the kicked Hamiltonian. 

Definition 6.9 (V-Smooth) Let V be a unitary operator. A is V-smooth if 
and only if for all f G Ti, V(t)cf) e D{A) for almost every t G M and for some 
constant C, 

< c ii^ip. 

n 

Theorem 6.10 If A is V-smooth, then R{A*) c 'HaciY)- 

Proof Since HadV) is closed, we need only show R{A*) c 'Hac{V). Let 
0 G D{A*), f = A*<p, and let dp^ be the spectral measure for V associated with 
ip. Define, for the period, T, in (6.1), 

RdT) = {A*d [V{T)rf ). (6.44) 

V 27r 

We calculate, dropping the T for clarity, 

V /TT 

< ^ 11.^11 piz-'V-ll. 

Because A is (/-smooth, we see that 

n n 


< oo. 
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Thus, e L^(R). By the Riesz-Fischer theorem (4.26 Fourier Series, p. 96- 
7,[100]), :F(e) e lK 

The speetral resolution of 1/[T] is 

p27V 

V[T]= e^UErie), 

Jo 

so 

I*27T 

{V[T])^= / E^UErie). 

Jo 

Therefore, from (6.44) we obtain 

1 

= ^ {A*(P,E^UET{e)ij) 

v27r Jo 
1 

VZTT Jo 

VZTT Jo 

Using the inverse of the expression above for JF(6*) gives 

En = / E^^Eie) d9. 

V27r Jo 

As we have just shown that E{9) G d^^{9) = E{9)d9 is absolutely eontinu- 
ous, whieh implies that 'll; e R{A*) is in Hac{V) and so R{A*) C Hac{V). □ 

Theorem 6.11 (Unitary Putnam-Kato theorem) Let V be a unitary 
operator, and A a self-adjoint bounded operator. If C = U[U*, A] > 0, then V 
is absolutely continuous on i?(C'^/^). Hence, if R{C^^‘^) is cyclic for V, then V is 
absolutely continuous. 

Proof. The diserete time evolution of an operator A is given by 


En = V-^AV^. 
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Calculate 


- Tn-x = C-MC" - !/-(«- 

= V-^[AV -VA]V^^-^'^ 

= V-^VIV-^A - AV-^]VV^-^ 
= V-^V[V*, AIV'^ 

= Gn, 


so 


G„4') = 5 ] (4>, V-"-V[V‘, A]V 

n=a n=a 

b 

n=a 
b 

n=a 
b 

n=a 

where C = V[V*, A]. We also have 

b 

n=a 

Taking the modulus and using the Schwartz inequality yields 

b 

y]||cW>||2<2|(,#,,l/-‘yll/V)| 

n=a 

= 2 I (CV, AV 

<2||A||||CV|P 

= 2 " /'iiii^ii2 
< oo 


and thus l/-smooth. 
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Finally, that V is absolutely continuous on R(C^^^) follows directly from The¬ 
orem 6.10. □ 


6.4 Finite rank perturbations 

Here, I utilise the results of Section 6.2 to show that perturbations of the form 
(6.5) lead to a pure point spectrum for the Floquet operator for a.e. perturbation 
strength A. 

I use directly the definition of strongly H-finite from Howland. 


Definition 6.12 (Strongly H-finite) Let H be a self-adjoint operator on 
7i with pure point spectrum, fin a complete orthonormal set of eigenvectors, and 
Hfin = Oinfin- A boundcd operator A ■. 7i ^ ICis strongly H-finite if and only if 




< oo. 


(6.45) 


n=l 


If H is thought of as a diagonal matrix on I2, i.e., H = \ fin){fin\, and 

A as an infinite matrix {aij}, i.e., A = n ^rnn \fim) {fin\, then (6.45) says 


E 


E 


< 00. 


(6.46) 


For our purposes, we need to show that if A is strongly iF-finite, then it is 
(/-finite. To satisfy the assumption that is trace class in Lemma 6.5 (and hence 
also compact in Theorem 6.4) we also need to show that A is trace class. 


Theorem 6.13 If A is strongly H-finite, then given U = for the period 

T in (6.1) and Hfin = O-nfiny 


(a) A is trace class, and 
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(b) A is U-finite. 

Proof, (a) Simply consider 

Tr(A) = ^{(t)i\A\(t)i) = ^aii< ^\aii\ 
i i i 

For each term in the sum (6.47) we trivially have 


\aii\ < 

and thus (6.47) is finite so A is trace class. 


(b) Noting that 




we calculate, by insertion of a complete set of states, 

n 

= _ u*e-^0-) (1 - 

= Y,{ct)n\A\ct)^){ct)^\ _ 

n m '■ ' ' 


A'\<l>n) 


E 

Ti,m 

E 

n,m 

E 

m 

E 




(1 - e-*^“ m 

_ {firn\A*\(t)n){(t>n\A\(f)rn) _ 

(1 - m _ fiiToi 

m /^e*®+) 

{firnlA* A\(j)rn) 

I g 6 giT" OLjfi / I 


I ^ g — E ^iT olji !Ti I 


2* 


The trace norm is then 


TrG.(9) = y; 


I^0ri 


^ _ g €. ^iT OLji j ft 


i2 * 


(6.47) 
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If this is bounded for e = 0, then it is trivially bounded for all e > 0. By (6.45) 
and a slightly modified version of (Theorem 3.1, [66]) this is finite a.e. for e = 0. 
Thus the trace norm of exists as e —0, which implies that the strong limit of 
Ge exists and we conclude that A is (7-finite. □ 


Theorem 6.14 Let U be a pure point unitary operator, and let Ai ,..., be 
strongly H-finite. Assume that the AkS commute with each other. Then for a.e. 
A = (Ai,..., A/v") i^ , 

1/(A) = e'‘{'^k=i^kAiAk)/hjj 

is pure point. 


Proof. This is a trivial modification of (Theorem 4.3, [66]). Let 


N 


/C = 0i((Alfc). 


k=l 


The elements of K, are represented as column vectors. The operator A : TC ^ IC 
is defined, for y eH,hy 


Ay = 

and therefore A* : IC His given by 


- 1 


Xi 

AnV 


Xn 


A*x = A\xi +-h A*^xisf. 

Accordingly, 0^(0) : /C —> /C, the matrix equivalent of equation (6.9), is intro¬ 
duced, 

G,i 9 ;U,A) 

= A[1 - [1 - Ue^^+] 
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The diagonal terms are finite a.e. beeause eaeh is [/-finite by Theorem 6.13. 
The off-diagonal terms are of the form X^X 2 , and so the Sehwartz inequality, 

ix*x2i'<iixiinix2i|2 

ensures that they are finite a.e. too. Henee, A is [/-finite as every term in the 
matrix G^{9; [/, A) is a.e. finite as e —> 0. 

The Hamiltonian may now be written as 

OO 

H{X) = Ho +A*W{X)Aj2^it-^T) (6.48) 

n=0 

and the Floquet operator as 

y(A) = 

where W (A) = diagjAfc}. In this form, the formalism of Section 6.2 is essentially 
fully regained, and we can proceed to apply Theorems 6.2, 6.4, 6.6 and 6.7. 

To establish the absolute continuity of the multiplication operator V on the 
space L^(R^; H) we proceed as in Proposition 6.8. Write 

V{G, ...An) = 


define 


where pk 


N 

D = — arctan(pfc/2), 

k=l 

—id/dtk, and compute 


N 

C = Y[Y*,D] = > 0. 

k=l 

In obtaining C as a direct sum of the Ck, we have had to assume that the AkS 
commute with each other. This complication comes when considering the term 


V{ti, ..., t]s[)V*(ti, ..., tk-i, Vkj tk+i, ..., /at) 
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in the equivalent of (6.34). To obtain the required form of require 

that the eommute.^ 

Moving on, each > 0 is equivalent to C in Proposition 6.8 and hence 
positive. Finally, we must show that is cyclic for V. This is no longer 

trivial as, for each k, while R{Cl^‘^) = R{A\), the range of A^. is not cyclic for U, 
hence V. To proceed, first note that 

k 

Now, as argued in Howland, we can assume that R{A*) is cyclic for U. To elabo¬ 
rate, define M{U, R{A*)) to be the smallest closed reducing subspace of H con¬ 
taining R{A*). If R{A*) is not cyclic for U, then HQ M is not empty. However, 
as shown below, ifyeHQ M, then A*W Ay = 0, so in ?-( 0 M., V{t) = U and 
is therefore pure point trivially. Thus, we can ignore the space HqA4 and restrict 
the discussion to M —i.e., we may assume R{A*) cyclic for U. 

The above relied upon showing that A*WAy = 0 for y E H Q AA. I now 
prove this. If y eHQ M and y' e M, then 

{y,y') = o. 

Given y' E AA, there exists anx E JC such that y' = A*x, so 

{y,A*x)=0. 


That is 


{Ay,x) = 0. 


This is true for all x E fC. Suppose y" E H. Then WAy" E /C and so 
_ {Ay, WAy") = 0. 

^This restriction is not required in Howland’s self-adjoint work because the summation over k 
in the Hamiltonian (6.48) enters directly, rather than in the exponent of V. 
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That is 

{A*WAy,y")=Q. 

As this is true for any y" G Ti, we conclude that A*WAy = 0 on Ti © Ad. 

Thus, R{A*) (with A acting on may be assumed cyclic for U, 

hence cyclic for V. 

I must finally show that = R{A*). We have 

/?(/!•) = |jfl(/ia = Ufi(c© 

k k 

and 

R{C) = [jR{Cu). 

k 

As R{A*) = R{A*A), R{C^/^) = R{C) and we have shown that R{C^/^) = 
as required. □ 

Finally, I wish to make the connection with my original aim—to show that 
Hamiltonians of the form 

N oo 

H{t) =Ho + J2 {'ipkl (6.49) 

k=l n=0 

have a pure point quasi-energy spectrum. 

Theorem 6.15 Let Hq be pure point, and define our time-dependent Hamilto¬ 
nian as in (6.49). Iffiii, . .., G Ii(Hq), then for a.e. A = (Ai, ..., X^) in 
the Floquet operator 

V = F(L^^=A>=^9k:)(^kl)/fifj 

has pure point spectrum. 

Proof. This theorem is just a special case of Theorem 6.14 with the A^s given 
by \'fk){'fk\- Noting (6.7), the A^s clearly commute. As Howland shows, |'0)(V’l 
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is strongly i7-finite if and only if '0 e Thus Theorem 6.14 applies and the 


result follows. 


□ 


6.5 Discussion of results and potential applications 

Of fundamental importanee in showing that the quasi-energy speetrum remains 
pure point for a.e. perturbation strength A was the faet that '0^ G Ii{Hq). That is, 
if we write 



n=0 

where the \(j)n) are the basis states of Hq, then ipk G h{Ho) if and only if 



If this requirement is dropped, and we only retain 'ip^ G then (Theo¬ 

rem 3.1, [66]) fails and there is the possibility that I^(A) will have a non-empty 
eontinuous speetrum. It was this faet that Milek and Seba [90] took advantage 
of in showing that the rank-1 kicked rotor could contain a singularly continuous 
spectral component under certain conditions on the ratio of the kicking frequency 
and the fundamental rotor frequency. They analysed two regimes of the perturba¬ 
tion. One where xp G in which case the numerical results clearly showed 

pure point recurrent behaviour, and the other where ^p G 12{Hq), but ip ^ Ii{Hq). 
In the second case, the authors further proved that the absolutely continuous part 
of the spectrum was empty,® and thus the system contained a singularly continu¬ 
ous spectral component. The numerical results reflected this, with a diffusive type 
energy growth being observed. 

turns out that Milek and Seba actually made an assumption in obtaining this result which is 
as yet is unjustified. Chapter 7 investigates this in detail. 
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With the generalisation of Combescure’s work here, namely Theorem 6.15, it 
is now possible to investigate the full class of rank-iV kicked Hamiltonians. A 
sufficient requirement for recurrent behaviour has been shown to be ^|Jk € 
and so I must turn my attention to perturbations where this requirement is no 
longer satisfied. This is the topic of Chapter 7. 

6.6 Summary 

I have shown, in a rigorous and general fashion, that the spectrum of the Floquet 
operator remains pure point for perturbations which are constructed from projec¬ 
tion operators that are in turn built from Hilbert space vectors which are elements 
of Ii{Hq). In simple terms, with the Hamiltonian perturbation W = \^p){^jJ\, one 
requires 

<00 

n 

for the Floquet operator to have a pure point spectrum for almost every perturba¬ 
tion strength. 

As alluded to in Section 6.5, to investigate systems that may display chaotic 
behaviour I would like to relax, in a controlled manner, the conditions that lead to 
a pure point spectrum for V. The emergence of a continuous spectrum is a vital 
ingredient in making further progress. Chapter 7 investigates this question. 
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7.5 Generalising the results of Milek and Seba. 128 

7.5.1 Discussion . 130 

7.6 Summary. 131 


This chapter extends the results of Combeseure [28] in a number of ways. 
Firstly, as an alternative to the work just presented in Chapter 6 I direetly extend 
Combeseure’s result on rank-1 perturbations to rank-iV perturbations—pleasingly, 
the same results as in Chapter 6 are obtained. I will then move on to the impor¬ 
tant question of the emergenee of a eontinuous speetral eomponent of the Floquet 
operator. I extend all of the results of Combeseure and also those of Milek and 
Seba [90] to rank-iV perturbations. 

The investigations lead, in a natural way, to a eonjeeture presented by Combes¬ 
eure [28] eonceming the dependence of her results on the partieular Hq eigen¬ 
value sequenee. The analysis herein leads to a number-theoretie eonjeeture that 
has stood for over fifty years [112]‘ on the estimation of finite exponential sums. 
Work already done in this area [15] will be examined in detail. 

In examining the work of Milek and Seba, I highlight a number of miseon- 
eeptions and reetify them. Worryingly though, their work is not in faet fully 
justified—a point so far missed in the literature. The resolution is direetly linked 
to the number-theoretie investigations just mentioned. 

'The reference is to the 1954 English translation of Vinogradov’s original work, published in 
1947. The work in Vinogradov’s 1947 monograph incorporates results from a series of papers 
and a first monograph from 1937. It is unknown (to me) when the conjecture I refer to was first 
presented, but it was at least fifty years ago. 









7.1 Outline and summary of results 


105 


7.1 Outline and summary of results 

In this chapter, I consider Hamiltonians of the form 

( N \ oo 

^5(t-nT) (7.1) 

fc=l / n=0 

where G M and each vector \^jJk) is a linear combination of the Hq basis states, 
\(j)n)- See Section 6.1 for the definitions and properties of these objects. 

As already discussed, the basic result is that if every \^jJk) G the spec¬ 

trum will remain pure point for almost every set of perturbation strengths A^. If 
this condition is dropped for any one of the |^fc) then Vai,...,Ajv is longer pure 
point. On the subspace Tik, the space for which |^fc) is a cyclic vector for the 
operator U, the spectrum is purely continuous. 

The other key result of this chapter concerns a number-theoretic conjecture 
stated by Vinogradov [112]. For Milek and Seba’s work to be properly justified, a 
sufficient condition is for Vinogradov’s conjecture to be true. This observation is 
linked to the conjecture put forward by Combescure [28] and partially addressed 
by Bourget [15]. 

In Section 7.2 I extend Combescure’s rank-1 theorem on the pure point spec¬ 
tral nature of V to the rank-A^ case. In Section 7.3 I then show the existence 
of a continuous spectrum for the case where 77o is the harmonic oscillator and 
the perturbation is rank-iV. In Section 7.4 I investigate Combescure’s conjec¬ 
ture, the partial answer provided by Bourget and the link to number theory and 
Vinogradov’s conjecture. Finally, in Section 7.5, 1 extend Milek and Seba’s work 
to the rank-A^ case. A number of conceptual and mathematical errors are firstly 
highlighted and then resolved. 
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7.2 A rank-N generalisation of Combescure’s first 
theorem 


Consider the measures 

mk,\^ = {^pk\Ex^{S)\^pk)^ 


Each l^fc) admits a cyclic subspace of H, Hk- As argued in the later part of the 
proof of Theorem 6.14, on the space H 0 Hk ^, the perturbation 

N 

fe=i 


is null and thus V = U is trivially pure point. Henceforth, we may safely restrict 
the proof to the subspace 0^=^ Hk for which the vectors 1^^) form a cyclic set. 

Directly following Combescure, the measure for a point x E [0, 2%) for the 
operator V acting on the state |'0fc) is given by 

'I Bk{x), (7.2) 

^^k 

where 


and 


Bk{x) = 


-27r 


dmk,\^=o{9) (sin^ [(x - 6')/2]) 


-1 


uo J 

This result is the essence of Lemma 1 in Combescure’s work. When Hq is pure 
point, it is a trivial calculation to show that 


B 


-i, 


X) = 


V- 

Cl 


I {flk)7 




n=0 


(7.3) 


Corollary 2 in Combescure’s work is replaced with the following. 
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Theorem 7.1 Assume Hq is pure point, with {(/>n}neN as eigen¬ 

states and eigenvalues. Let each 


n =0 

be cyclic for Hq (hence, cyclic for U and V) on Tik and ('0^1'^^) 
belongs to the point spectrum o/VAi,,,,,AAr if and only if 

N 

fc=i 

where 

On = 2'K{an/‘l'Kh}, 


Ski- Then e“ 


{z} being the fractional part of z. 


Proof (7.1) The proof follows that in Combescure. By the cyclicity of each 
on Tik and the argument in Theorem 6.14, e“ is an eigenvalue of Vai,...,A]v if 
and only if every rnk^\^{{6}) 7 ^ 0 at 0 = x. As already mentioned, using 


d'rnk^\f.=Q ^ ^ 1 ( 0^)711 d(d OjfjdO 


n=0 


we obtain, for each k. 


Bk\x) = 


I i^^k)r 


^ ^iv?[{x-en)/2\ 


n=0 


Now consider the eigenvalue If it were to be that for some k, rnk,\^. ({a;}) = 0, 
then we would have found a vector, namely \fk), such that V\^jJk) was continuous. 
We have in fact found that the whole subspace TLk is continuous. Thus, for V to be 
pure point, every mk,\f. ({ 6 *}) 7 ^ 0. Thus, we are lead to consider the requirement 

N 

< 00. 

k=l 


□ 
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As in Combescure, the relationship 

^ |(afc)npcotg = cotg^ (7.4) 

also holds for each k. To show (7.4), consider each k separately. The proof is the 
same as for the rank-1 case. See [28]. Points to consider are that each projection 
operator in the rank-TV projection is normalised and hence for every k we have 

N 

n=0 

In order to complete the generalisation of Combescure’s first theorem, we require, 
just as in Combescure, two additional Lemmas. 

Lemma 7.2 IfYl^=o l(®fc)n| < C) 0 , then < oo for almost every x G M. 

For each fc G 1,..., the proof is identical to that in Combescure. 

Lemma 7.3 The following two statements are equivalent. 

(a) For almost every (Ai,..., Aat), Vai,.,.,Ajv only a point spectrum. 

(b) For every k E {1,..., N} and for almost every x, Bk{x) 0. 

The proof is again virtually identical to Combescure’s proof. For each k, the 
continuous part of the spectrum is supported outside the set = {x E [0, 27r) : 
Bk{x) 7 ^ 0} and, for A^ 7 ^ 0, the point part of druk^x,^ is supported by the set Ek. 
Thus, for to be pure point for almost every Ai,..., Aat and for every k, 

we require 

mfc,Afc([0,27r) \ Ek) = 0. 

This in turn implies that for every k 

2tt 

dX'k /i(Afc)mfc,Afc([0, 27r) \ Ek) = 0, 
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where = Xk/h and 


h{X) = 23fJ 


1 


1 — ce*^ 


for some |c| < 1. 

Combeseure’s Lemma 5 trivially applies for eaeh k. Thus, I have generalised 
Combescure’s work to obtain the result that the Floquet operator for the rank- 
N perturbed Hamiltonian has a pure point spectrum. As already mentioned, the 
result matches that obtained in Chapter 6. 


7.3 A rank-N generalisation of Combescure’s sec¬ 
ond theorem 

Having shown that the Floquet operator remains pure point for perturbations con¬ 
structed from the vectors \^jJk) € Combescure relaxes this condition to 

allow for the emergence of a continuous spectral component of the Floquet op¬ 
erator. This result is easily generalised to the rank-TV case. The key point is that 
Combescure’s technique applies independently for each k. I do not discuss the de¬ 
tails of the rank-1 proof here at all, delaying an analysis to Section 7.4 where I will 
have the opportunity to generalise the results still further. Here, I simply provide 
the argument for why each k may be treated independently. Before proceeding, 
some subtleties of what Combescure actually shows are highlighted. They are 
seemingly overlooked by some in the literature (e.g., [90]). 

The cyclicity requirement was essential in the proof that the Floquet operator 
spectrum was pure point. Here, we can happily ignore the cyclicity conditions, 
as our only goal is to establish the existence of a state in the continuous subspace 
kiconi- We need not try and ensure the result obtained by considering (^|i7(5') \'ip) 
is applicable to all other vectors in 77—the very idea is ill-formed as the pertur¬ 
bation is null on a subset of 77 and thus there is always part of 77 where V has 
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a discrete speetrum. Milek and Seba seem to have missed this point, restating 
Combeseure’s theorem in a way that implies that all '0 are in TYcom- 
If 

{il)k\E{{x})\ipk) = 0 

then Tiac contains at least the state \ipk)- The point to be mindful of is that this does 
not allow one to conclude that the Hilbert space for the operator V has Hpp = 0. It 
is this error that Milek and Seba [90] have made. To draw that eonelusion would 
require an argument to show that a eyclic veetor does in fact exist for V. This does 
not seem possible in the general eontext here. 

Combeseure’s proof (Lemma 6 in [28]) that crcont(C) 7 ^ 0 is based on showing 
that B~^(x) —>■ 00 (equation (7.3)). As the speetral measure of a single point x is 
proportional to B(x) (equation (7.2)), if B~^(x) —> cx), then the eontribution of 
the single point is zero. That is, is in the eontinuous spectrum of the Floquet 
operator. Combescure argues (see Section 7.4 for details) that 

B~^(x) > #S'(x) 

where ^S{x) is the number of elements of a partieular set S. She then shows (the 
bulk of the proof) that 4 ^S{x) 00 and thus B~^{x) 00. I generalise the 

result in a straightforward manner. 

Theorem 7.4 Assume = nhw with lo irrational. If I'fk) ^ h{Ho) for at 
least one k E 1,... ,N, then acontiV) 7 ^ 0. 

Proof. (7.4) Following the same argument as for the rank-1 ease, we take 

|(afc)n| = rr^<2'K 

for the state \iJk), in such a way that the eondition {f)k\'fi) = ^ki is preserved. 

With this eonstruetion, Combeseure’s proof that the number of elements in 
S{x) is infinite applies to eaeh subsequenee Sk{x). The number of elements. 
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in each sub-sequence for which l'^^) ^ is infinite. The Floquet 

operator for the rank-iV perturbed harmonic oscillator obtains a continuous spec¬ 
tral component. □ 


7.3.1 Discussion 

It must be noted that the proof presented by Combescure (Lemma 6, [28]) is only 
valid for the eigenvalue spectrum, 


On = nhuj 

of the harmonic oscillator. Combescure does however conjecture that the argu¬ 
ment will be valid for more general eigenvalue spectra, including the rotor, 

an oc 

For Milek and Seba’s numerical work (using the rotor) to be based on valid math¬ 
ematical arguments, a proof of this conjecture is required. Currently, no such 
proof exists. In Section 7.4 I show that if a conjecture from number theory on 
the estimation of exponential sums is true, then Milek and Seba’s work can be 
justified. The rank-iV generalisation is straightforward. Considering the number 
theory conjecture has stood for over fifty years, it seems we may have to wait quite 
some time for a proof. 

For more general eigenvalue spectra (loosely oc n^) the situation is signifi¬ 
cantly better. For J > 3 Bourget [15] has made significant progress. A continuous 
component of the Floquet operator exists for certain constructions of |V’)- The 
conditions are complicated and more restrictive than the |'0) ^ h{Ho) condition 
for the harmonic oscillator. The result is easily extended to the rank-iV case due 
to the independence of each k as already discussed. Utilising the same number- 
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theoretic conjecture as in the j = 2 case will also allow for improvements to the 
work of Bourget. See Section 7.4. 

Returning to the harmonic oscillator case, by applying Theorem 7.4 we may 
conclude that for each 1^^.) ^ Ii{Hq), Hk is purely continuous. Thus, by dropping 
the /i condition for all I have shown that V is purely continuous on the 
subspace of H where the perturbation is non-zero. On the subspace of H where 
the perturbation is zero, V = U trivially and thus that portion of the Hilbert space 
remains pure point. 

7.4 Combescure’s conjecture and number theory 

Combescure makes a remark (Remark c., [28]) that she believes Theorem 7.4 
(Lemma 6, Combescure [28]) is generalisable to include systems other than the 
harmonic oscillator. Explicitly, she conjectures that Hamiltonians, Hq, with eigen¬ 
values, an, of the form 

V 

an = h'^(3jn? (7.5) 

i=o 

with (3jT j2n Diophantine for some j : 1 < j < p will have the vector ^ in the 
continuous spectral subspace of 14 . 

At an intuitive level, one would expect this to be true. The precise nature of 
the eigenvalue spectra (proportional to n or a polynomial in n) should not make 
a significant difference. Milek [89] argues that Combescure’s work can be used 
in the case based on evidence from some numerical work that shows that the 
sequences obtained are “almost random”—however, the argument is not entirely 
convincing to me. The cited numerical work of Casati et. al. [20] discusses the 
existence of correlations in the energy levels, rather than the lack of correlations. 
While the deviations from a Poisson distribution look small to the naked eye, 
Casati et. al. [20] find deviations from the expected Poisson distribution of up to 
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17 standard deviations. The energy levels are correlated—it is arguable that they 
are not characterisable as “almost random” as Milek asserts. 

I began to explore the possibility of developing a proof to the conjecture. A 
few interesting results have come from this investigation and will be presented 
here. While doing this work, I was unaware that in late 2002, Bourget [15] pro¬ 
duced a proof of a slightly modified conjecture for all but the p = 2 case in (7.5). 
The techniques used by Bourget are the same as those followed in my work. I will 
analyse Bourget’s work, and highlight the key breakthrough made. I also provide 
a modified argument to obtain the proof which is, I believe, significantly easier 
to follow. Importantly, it also covers the p = 2 case missed by Bourget due to 
technical difficulties. However, it comes at the expense of relying upon a (quite 
reasonable) conjecture. I do not claim that what is presented is adequate on its 
own, but it does play a complementary role in understanding, or perhaps appre¬ 
ciating, Bourget’s proof. The reliance on the conjecture simply removes the need 
for much of the technical wizardry in Bourget’s proof. Use of the conjecture also 
strengthens the work. The work also indicates, or highlights, that Combescure’s 
conjecture is solved by a number-theoretic conjecture that has stood for over fifty 
years. What seems a perfectly reasonable conjecture on physical grounds is shown 
to be directly related to an abstract mathematical conjecture. 

In what follows, I will rely heavily upon the lemmas and theorems in Chapter 2 
of [84]. I also use some results on Weyl sums from [112]. Of key importance is 
an understanding of Combescure’s proof of her Lemma 6 on the emergence of a 
continuous spectrum for the kicked harmonic oscillator. This will be discussed at 
the appropriate time in this section. 
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7.4.1 Number theory 

To investigate Combeseure’s eonjeeture we require two eoneepts from number 
theory—the elassifieation of irrational numbers and the discrepancy of a sequenee. 
I first introduee the concepts and define the relevant ideas. I then proceed to anal¬ 
yse the conjecture and the proof provided by Bourget. As the discussion pro¬ 
gresses, the new work that I have done will be presented. 

For any number (3, define 

• [/3], the integer part of (3, 

• {(3}, the fractional part of (3, and 
. (/3)=min({/3},l-{/?}). 

{(3) is simply the “distance to the nearest integer”. Definition 7.5 is taken directly 
from Kuipers and Niederreiter (Definition 3.4, p. 121, [84]). 

Definition 7.5 Let rj be a positive real number or infinity. The irrational, (3, is 
of type p if p is the supremum of all r for which 

lim q^{q(3) = 0, (7.6) 

n—^■oo 

where q runs through the positive integers. 

The idea behind this definition can be seen by considering rational f3 = p/q' 
for integers p and q'. Run through the positive integers q. At q = q', (qfi) = 0, and 
so there is no supremum p for r in (7.6). In effect, p ^ oc. For irrational (3, (qP) 
is never equal to zero but will approach zero. If the approach is very slow, then 
a small r is enough to prevent (7.6) from approaching zero. {qP) approaching 
zero slowly is, in a sense, indicative of P being badly approximated by rational 
numbers. Even for very large q', p/q' remains a poor approximation to p. Thus, 
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the smaller r], the stronger the irrationality of (3. This is reasonable in the sense 
that rational (3s act like numbers with 77 —> cxo. As stated in [84], all numbers (3 
have type r/ > 1. 

I now define the discrepancy of a sequence—a measure of the non-uniformity 
of the sequence. Consider a sequence of numbers^ Xn in [0,1), 

OJ = (a:n)neN with G [0,1). 

For 0 < a < 6 < 1 and positive integer N, A{[a, b), N) counts the number of 
terms of the sequence (up to xn) contained in the interval [a, b), 

A{[a, b),N) = #{n < N : Xn ^ [a, b)}. 


Definition 7.6 The discrepancy of the sequence uj is 

A([a,fe),iV) 


Dn{u) = sup 

0<a<&<l 


N 


-{b-a) 


(7.7) 


If the sequence 00 is uniformly distributed in [0,1) then ^ 0 as N —>cx). 
In this case, every interval [a, b] in [0,1) gets its “fair share” of terms from the 
sequence u. 

Estimating the discrepancy of a sequence will turn out to be vital in the anal¬ 
ysis of Combescure’s work. The sequence of interest is basically the eigenvalue 
sequence for Hq, but I will discuss this in greater detail later. 

The starting point for the estimations that we require is (equation (2.42), Chap¬ 
ter 2, [84]). This is a famous result obtained by Erdos and Turan. It states that 


Dn<C 


m 


E 

h=l 


h 


N 


N 

E' 

n=l 


^2'KihXn 


(7.8) 


^Equivalently, consider any sequence Xn and consider the discrepancy of the sequence mod¬ 


ulo 1. 
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for any real numbers xi,..., xat and any positive integer m. The sum 

N 

^2TTihXn 

71=1 

is an example of a elass of exponential sums known as Weyl sums, refleeting the 
pioneering work of Weyl on providing estimations for them. Vinogradov [112] 
improved on some of the estimations of Weyl. Weyl and Vinogradov’s results 
eoneern the modulus of the sum, \S\, and characterise it as 

\S\<7N, 


where N is the number of terms in the sum and 7 tends to zero as iV —> cx). The 
subtle behaviour of 7 is linked to the rational/irrational nature of the terms in the 
sequence. 

I will use some basic results from the introductory chapter of [112]. In general, 
write 

N 

S = exp {2mF{n)) 

n=l 

for some function F{n). The application here is when 

F{n) = f3n^. 

For (3 rational {not the case I will be interested in) L. K. Hau proved that IS"! was 
of order 

7Yi-(i/i)+e 

(p. 3, [112]) and that this estimate could not be much improved. Here, I am 
interested in the case where /3 is irrational. Estimations are much more difficult, 
and form the major aspect of the work by Vinogradov. The estimations depend 
upon making a rational approximation to (3 and are complicated functions of N 
and j. Very loosely, he obtains results like 


1^1 = 0{N^-P') 
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where 

^ 3(j - 1)2 log 12j{j - 1) ■ 

Vinogradov states 

It is a plausible eonjeeture that the estimate in (7.9) holds with p' 
replaeed be 1/j — e ... A proof or disproof of this conjecture would 
be very desirable. 

With the dual aim of extending Bourget’s proof to the p = 2 case (the rotor consid¬ 
ered by Milek and Seba) and “simplifying” Bourget’s proof, I state this conjecture 
formally. 

Conjecture 7.7 Consider the sum 

N 

S = exp 27 rin^ / 3 j. 

n=l 

For all N greater than some critical value, 

1^1 < 

for all e > 0 and some constant c G M. 

I do not attempt to prove Conjecture 7.7. Given the lengths gone to by Vino¬ 
gradov to obtain the results presented above, it seems rather unlikely that a proof 
or disproof will be found any time soon.^ 

^Incremental improvements on the estimations presented by Vinogradov in [112] have been 
made over time. While Bourget [15] makes use of these improved results, the conjecture itself 
remains unproven which is the only result of any consequence in this discussion. 
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7.4.2 New results on discrepancy—upper and lower bounds 

Armed with the estimations on Weyl sums, I now proeeed to derive both upper 
and lower bounds on the discrepaney for sequences of the type 


(n^/3) 


for [3 of any type r; > 1. It must be remembered that the upper bound obtained 
is contingent upon Conjecture 7.7. The lower bound obtained is not dependent 
upon any unproved conjectures. The result obtained highlights the “best possible” 
nature of the conjectured upper bound. 

Firstly, (Lemma 3.2, p. 122, [84]) is generalised to arbitrary j. 

Lemma 7.8 The discrepancy D^iujj) ofooj = {n^(5) satisfies 



for any positive integer m and e > 0, where C and c' are absolute constants. 

Proof. (7.8) Consider equation (7.8). It is applicable to the first N terms of the 
sequence ooj. We have 



(7.10) 


for any positive integer m. Consider the sum over n. 



Conjecture 7.7 allows this sum to be bounded by 




We are free to write 


sin7r/i/9| 


c = 
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as sin nh(3 is just some positive real number. Substituting this result into (7.10) 
gives 



Now following the argument at the end of (Lemma 3.2, [84]) the desired result is 


obtained. 


□ 


I now give the generalisation of (Theorem 3.2, [84]). It provides the “best” 
upper bound one eould hope for when estimating the discrepancy of the sequence 
ojj = (n^/3). Again, remember that the proof relies on Conjecture 7.7. 

Theorem 7.9 Assume Conjecture 7.7 is true. Let jS be of finite type r]. Let j 
be a positive integer j > 1. Then, for every e > 0, the discrepancy Djq^ujj) of 
ujj = {n^(3) satisfies 


DN{ujf = O . 


Proof. (7.9) Let e > 0 be fixed. By (Lemma 3.1 and Lemma 3.3, p. 121- 


3, [84]), 



for a fixed e' > 0. Combining this with Lemma 7.8 gives 



for all m > 1. Now choose m = j. We obtain 

DN{ujj) < C 


where e = e" + e'/pj. 


□ 
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Theorem 7.9 is, in a sense, optimal. For functions f,g, define / = fl{g) if 

f/g ^ 0 . 

Theorem 7.10 Let (3 be of finite type r]. Let j be a positive integer j > 1. Then, 
for every e > 0 , the discrepancy D^iuj) ofojj = {n^(3) satisfies 

Dn{uJj) = n . 

Proof (7.10) Let e > 0 be fixed. For any given e' > 0, there exists 0 < 
6 < p with l/{p — 6) = (l/p) + d. By (Definition 3.4, p. 121, [84]) we have 
hmg^oo {q/3 — j) = 0 and thus 

{qp) < 

for an infinite number of positive integers q. There are infinitely many positive 
integers q and p such that 

\P-Plq\ < 

That is, by choosing q large enough, we can always find a p such that \qP — p\ = 
(g/3). As q increases p/q is a better approximation to the irrational fi. For 6 some 
irrational with | 0 | < 1, we have 

/3 = p/g + 0 g-l-^+W 2 ). 

Pick a q such that the above relations are valid. Set 

N = . 

Then for 1 < 


n^(3 = n^{p/q) + 6^, 
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with 


^ ^-l-(< 5 / 2 )_ 


Thus, none of the fractional parts {/3}, {2^f3}, ..., { /3} lie in the interval 

J = so 


D^{ujj) > 


A{J,N) 

N 


A(J) 


A(J), 


where A( J) is simply the “size” of the set J. For large enough q, we have A( J) > 
l/2g. But from the definition of N it is clear that 


N < qA^-^) <N+1<2N, 


so 

q-^ > cNAAri-S)]-\ 

Combining these inequalities, we obtain 

Dpfiujj) > c'nAAv-S)]-^ 

where e = e'/j. That is, we have shown, for all e > 0, that 

DN{ujj) = n . 

□ 
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7.4.3 Combescure’s conjecture, Bourget’s work and new re¬ 
sults 


Before diseussing the eonjeeture, we must elearly understand Combeseure’s proof 
for the harmonie oseillator ease. As stated in Seetion 7.3, the aim is to show that 


n=0 


sin (x — 9r, 


oo. 


Define the set S'(x), 


S{x) = {n : \x — 9n\ < |an| = n '^27r}. 


(7.11) 


Each n is an element of S{x) if x is “close to 6*„”. Note that 9n = 2 ti{ an 
where {.} is the fractional part, not “set” and are the eigenvalues of the base 
Hamiltonian Hq. 

Given that sin a; < x for all x > 0, Combescure obtains a lower bound for 

2 E a •»w- nra 

n£S(x) 


Each n G S'(x) gives a contribution to the sum of greater than one as |a„|/|a: — 
9n\ > 1. By only considering ^S{x), we simply count a “1” each time. 

The results on the discrepancy of sequences are now used, with the sequence 
a.'HO = {9n/2'n). Note that each element of the sequence cuho is in [0,1). 

Consider the interval, defined for every x G (0, 27r) and centred around x / 2 t [. 




(7.13) 


Eor large enough N, Jn{x) C [0,1). The size of the interval is 2N . Using 
this particular subset and noting that the definition of discrepancy (7.7) involves 
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taking the supremum over all subsets of [0,1), Combeseure obtains 
\N-^A{Mx),N)-2N-^\ < D^iuuo). 

Multiplying through by N gives 

\A{Mx),N) - 2N^-'^\ < ND^ioono). (7.14) 

As lij) i hiHo) 

^ ^ |®n| ^ 

and thus 

1/2 < 7 < 1 

from simple eonvergenee arguments. Therefore, N^~'^ grows at a rate"^ less than 
At this stage, Combeseure utilises the theorems diseussed above on the 
diserepancy of sequences. For the eigenvalue sequence, a„ = nhu, of the har¬ 
monic oscillator^ the j = 1 case of Theorem 7.9 applies which is exactly (Theo¬ 
rem 3.2, [84]). Combeseure obtains the result^ 

D^{u;no)=0{N-^/^+^). 

For the sequence cuho, /? = a;/27r. If f3 is an irrational of constant type (rj = 1), 
the strongest type of irrational, then 

_ I^n{ooho) = 0{N^). 

"^Interestingly, it can in fact not grow at all (7 = 1) which is a subtle point seemingly missed 
by Combeseure and others. The rank-1 projection operator from the vector jr/;) constructed with 
7 = 1 is not shown to lead to the emergence of a continuous spectrum. Therefore, the statement 
that \tj}) ^ li{Ho) implies \ijj) € Ticont is not in fact proved to be true. There are states not in 
li{Ho) that may not be in the continuous spectrum. In practice (numerical, experimental work) 

this should not cause any trouble. It is clearly easy to avoid 7 = 1 . 

^Do not confuse w, the harmonic oscillator frequency, with ujho, the label for the sequence in 

[ 0 , 1 ), the discrepancy of which is being bounded. 

®This is not based on a conjecture as for j = 1 a direct proof is possible, bypassing Conjec¬ 


ture 7.7. See [84]. 
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As the right-hand side of (7.14) ean be made to grow arbitrarily slowly, we eon- 
elude that the left-hand side must grow slowly too. Thus, to eaneel the growth of 
2N^~'^, A{Jn{x), N) must grow at a rate arbitrarily elose to that of We 

see that 

A{J]^{x), iV) —> cx) 

as > cx). It is now a simple observation [28] that this implies that #*S'(x) —> cx 
and thus B~^{x) —> cx. Thus, is in the eontinuous spectral subspace of the 
Floquet operator V. 

The importance of the eigenvalue sequence is seen in that if we cannot limit 
the right-hand side of (7.14), then we cannot place a lower limit on A{Jn{x), N) 
and thus we cannot conclude that B~^{x) —*■ oo. Two barriers to limiting the 
right-hand side of this equation exist— j and r]. If, still in the harmonic oscillator 
case, we wished for P = ui/27 t to only be of a weaker type, say 77 = 2 , we would 
no longer be able to conclude that B~^ 00 . The right-hand side would grow 

like which is always faster than for 1/2 <7 < 1 which grows 

at a rate of Thus, no suitable lower limit for A{Jn{x), N) can be found. 

Similarly, if the eigenvalue sequence is generalised (Combescure’s conjecture) 
then we run into trouble. For j = 2, the lowest possible growth rate for the right- 
hand side we can obtain, taking Conjecture 7.7 as true, applying Theorem 7.9 and 
noting Theorem 7.10 which says we cannot do any better, is, once again, 

For larger j, the situation only gets worse. 

Given these seemingly significant problems, the natural question to ask is: 
“How does one get around this problem?”. The answer is provided in the work 
of Bourget [15]. Bourget proves a weaker theorem than Combescure’s conjecture. 
Where Combescure kept the same requirement on It/), that it be in Bourget 

has a j-dependent requirement. Essentially, for increasing j the an terms used to 
construct \ip) must decrease more slowly with n. See Bourget’s work [15] for the 
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exact requirement, which depends on the best estimates available for Weyl sums 
discussed earlier and thus is a non-trivial function of j. 

The key insight in obtaining the proof is to modify the set S{x) (equation 
(7.11)) and the corresponding interval Jn{x) (equation (7.13)) that are considered. 
Importantly, they become j-dependent. Bourget reduces the shrinking rate of the 
set Jn{x) as a function of N just enough so as to allow the weaker limits on the 
discrepancy to be good enough to force the right-hand side of the equivalent to 
(7.14) to be less than the left-hand side, while keeping strong enough control on 
terms in the new set S{x) to still argue that B~^ —> oo. 

Using the best available estimations on Weyl sums and plugging these into the 
upper bound formulas for discrepancy (as discussed earlier when introducing the 
work by Vinogradov), Bourget manages to provide a rigorous proof of the exis¬ 
tence of a continuous spectral component of the Floquet operator (the essence of 
Combescure’s conjecture) for j > 3, leaving only the j = 2 case unresolved. The 
proof is, unfortunately, unavoidably clouded by the “messy” estimates available 
for Weyl sums and thus, the essence of the proof is difficult to see. Here, I will re¬ 
visit the proof, but (utilising Conjecture 7.7) apply Theorem 7.9 which says (using 
2e, rather than e for technical reasons), for all e > 0 

With this very clean estimate, it is far easier to see how Bourget’s work provides a 
proof that a continuous spectral component of the Floquet operator exists. It also 
extends the result to j = 2. Of course, the j = 2 case remains unproved as I have 
relied upon Conjecture 7.7, but I highlight the fact that a solution to Vinogradov’s 
conjecture would solve Combescure’s physics conjecture. I have also simplified 
the j-dependence of the a„s used to construct \'ip). 


Theorem 7.11 Assume Conjecture 7.7 is true and thus Theorem 7.9 follows. 
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Assume (3 is irrational and of type p. Then for all positive integers, j, the Floquet 
operator, V, has (JcontiV) 7^ 0 < 7 < 1/2 + l/2pj. 

Proof (7.11) The proof relies upon the techniques utilised by Bourget. In 
essence, we simply increases the size of the interval (equation (7.13)) from 2N~'^ 
to Qog . The important change is the first factor. The log N term 

is essential for technical reasons, but has a negligibly small effect on the shrinkage 
rate of the interval for large N. As \ogN/N^^ —»■ 0 as iV —*■ cx) for all <5 > 0, for 
N large enough we have 

2JY2(l/2-7)QQgjY^-l/2 ^ 2iV2(V2-7-5)_ 

Using this underestimate for the size of the interval, we easily obtain the equiva¬ 
lent of (7.14), 

|A(Jjv(x),iV) < NDn{uJj), 

for the sequence tUj = (3). Now, using Theorem 7.9, it is evident that to ensure 

A{Jm{x), N) —>• cx), we must have 

2(1 - 7 - 5) > 1 - {l/pj) + 2e, 

or 

7 < 1/2 + (l/ 2 ? 7 j) -e-6. 

The condition 

1/2 < 7 < 1/2 + (1/27J), (V.15) 

where the “<” sign has absorbed the arbitrarily small numbers e and 5, must be 
satisfied to force A{Jn{x), N) —> cx). 

Finally, we must show that B~^{x) —> cxo when this larger interval is used. 
Corresponding to the new interval Jn{x), we introduce the new set S{x), 


S{x) = {n : \x — 6n\ < 27riV^^^'^^ logiV . 
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The estimate (7.12) is the same, exeept with the new set S{x), whieh no longer 
has all terms greater than unity. Thus, it is not enough to simply eount the number 
of terms in S{x). A more subtle estimate is required. Replaeing the numerator, 
|a„|, with something smaller, N~'^, and the denominator, {x — 9n), with something 
larger, logA^“^/^, we obtain 



n£S(x) 


which is essentially the estimate Bourget obtains. The estimate contained therein 
(Lemma 3.5 in [15]) then shows that B~^{x) —> oo and the argument is complete. 
□ 

Examining (7.15) note that for j = I (for rj = 1) we recover the simple 
result of Combescure. For all j > 2, we have a stronger (j-dependent) condition 
on 1^) than simply |^) ^ This complication is the main weakening of 

Combescure’s conjecture that Bourget and I have been forced to make. Note that 
the restriction on 7 takes into account the end point subtleties referred to in the 
preceding discussions. 

I have replaced the requirement that 1-0) ^ h{Ho) (i.e., 1/2 < 7 < 1) with 
the j-dependent requirement 1/2 < 7 < l/2 + (l/2j). In Bourget’s work, the re¬ 
quirement is stronger—directly related to the replacement of the known limits on 
Weyl sums (in terms of p in the earlier sections) with the “best possible” estimate 
from our Conjecture 7.7 of (1 /j) — e. 

7.4.4 Summary 

Reliance on Conjecture 7.7 and the result of Theorem 7.9 derived from it has 
allowed me to discuss Bourget’s proof without the complications of the messy 
estimations on Weyl sums. This simplified discussion highlights the key aspects 
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of Bourget’s proof. It has also shown that the j = 2 case for the emergence of a 
continuous spectral component of the Floquet operator is solved by Vinogradov’s 
conjecture. A proof of Vinogradov’s conjecture is no longer just of mathematical 
interest. It has a direct mathematical physics consequence. 

Finally, note that the rank-iV equivalent of this work follows in the same way 
as presented for the harmonic oscillator case in Section 7.3, providing a complete 
rank-7V generalisation of the work of Combescure [28]. 

7.5 Generalising the results of Milek and Seba 

Having established that the continuous subspace of H, Hconu is not empty, I now 
wish to characterise it—by identifying the singular and absolutely continuous 
components. Here, I extend the result of Milek and Seba to rank-iV perturba¬ 
tions. I do not extend the numerical results of Milek and Seba as they rely on 
the assumption that the j = 2 eigenvalue spectra lead to a continuous Floquet 
spectrum—a result I have just shown to be as yet unjustified. 

Theorem 7.12 Assume H{t) is given by (7.1) and that (7.2) applies. Assume 
B^^(x) — > oo and thus 'Hcont 7 ^ 0 - Then TCac = 0 and thus disc is not-empty. The 
Floquet operator, V, has a non-empty singular continuous spectrum. 

Proof. (7.12) As shown in the proof of Theorem 6.4a and easily calculated, 
the Floquet operator can be written in the form 

N 

V = + 

k=l 

where 

Rk = {e^^^'^-l)\fu){fk\U. (7.16) 

We can now use either (Theorem 5, Howland [65]) or (Theorem 1, Birman 
and Krein [8]). The theorem from the paper of Birman and Krein is more direct. 
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so we use it here. It states that if we have two unitary operators, U and V, that 
differ by a traee elass operator, then the wave operators 

= s- lim V^U-^P^,{U) 

I'^zLoo 

exist and their range is the absolutely eontinuous subspace of V, 

R{n±) = n.c{v). ( 7 . 17 ) 

We must show that the difference V — U is finite. With the notation in Chapter 6, 
where the perturbation W is given by ^4*^4 and 

A = \ip){'ip\. 


with 


we obtain 


1-0) = 


AiA*A = AiA = Y,{MM^i) 

I 

l,m,n 

l,m,n 

I 


as 1^) G hi^Ho) and is normalised. The perturbation to the Hamiltonian is trace 
class. The difference in unitary operators, U and V, is also trace class. By the 
triangle inequality for norms. 
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As ||f/||t,= 1, 

Tr (Eflt) < E 111 

\k=l / k l,m,n 



k 


= ^ a / 2(1 - cosXk/h). 

k 

Armed with a trace-class perturbation, we conclude that the wave operators exist. 
The existence of the operators il± means that they are defined for all states in the 
Hilbert Space H. Note (equation (7.17)) that the subspace 7fac(^) is equal to the 
range of these operators. However, Pac{U) gives zero when acting on any state in 
H because U is pure point. Thus, HadV) is empty. As 7-(cont is not empty, Hsc 
must be non-empty, and we have proved that a singular continuous subspace of 
the Floquet operator V exists. □ 

The key assumption in Theorem 7.12 is that > oo. This is certainly 

true for j = 1 if \ipk) 7^ For j > 2 the results were discussed in detail 

in Section 7.4. For j > 3, Bourget showed that one can construct vectors \dk) 
for which i?^^(x) — cxd. I have shown, in Conjecture 7.11, that if Conjecture 7.7 
is true then this result extends to j >2 and with improved requirements on the 
states Idk)- 

7.5.1 Discussion 

Milek and Seba make a number of incorrect statements in obtaining this result for 
the rank-1 case. Firstly, they state that the operator^ i? = [exp{i\/h) — 1] \d){d\U 
is rank-1 which it is not—the presence of the unitary operator U stops R from be¬ 
ing rank-1. This is not, however, important. The applicability of the theorems 
^As we are dealing with the rank-1 case, the subscript k may be dropped from (7.16). 
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in [8,65] does not rely upon the rank of the operator R, but upon it being of trace- 
elass. Seeondly, they elaim that the existence of the wave operators implies that 

a^cV (Z (7.18) 

This is, again, not true. Given that is empty, it is indeed possible to con¬ 

clude that (Jac(l^) is empty, as discussed above, but the relation (7.18) does not 
follow. Consider the situation where crcont((7) is not empty. Then there is a set of 
vectors in R. which are continuous for U. These vectors form the domain for the 
operator in the wave operators. The action with V does not however keep 
us in the subspace HconiiU) as the space we get to (the range for V^) is only in¬ 
variant for HcontiV), not HcontiU). Thus, we may obtain a vector, necessarily in 
7fcont(^) due to invariance, but possibly in Hs{U), and thus, we cannot conclude 
that crac(17) c <Jac{U). These two points discussed do not make the final results of 
Milek and Seba wrong, but “only” the proofs. 

Of greatest concern is the use of (Lemma 6, Combescure [28]) without jus¬ 
tification. Milek and Seba have assumed that Combescure’s conjecture is true. 
Bourget’s demonstration that a continuous spectral component of the Floquet op¬ 
erator exists does not cover the j = 2 case which is exactly the situation in Milek 
and Seba’s paper. Furthermore, I have shown that, using the “best possible” cases 
for discrepancy, the j = 2 case is covered, but, as I relied on Conjecture 7.7, I 
have not actually proved it. It is worrying that Milek and Seba’s work remains 
unjustified. 

7.6 Summary 

I have generalised the work of both Combescure [28] and Milek and Seba [90] 
from rank-1 to rank-A^. I have also discussed in detail Combescure’s conjecture. 
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my work on estimations of discrepancy and the demonstration by Bourget [15] 
that a continuous spectral component of the Floquet operator does exist for certain 
constructions of \ip). This covers the essential aim of Combescure’s conjecture on 
the existence of a continuous spectral component. A clear view of the essence 
of Bourget’s proof has been provided by taking a reasonable number-theoretic 
conjecture to be true. With this clear view, the work of Bourget becomes more 
accessible. I also demonstrated that reliance on Vinogradov’s conjecture allows 
one to extend the work to the j = 2 case, showing that a proof of Vinogradov’s 
conjecture would have direct implications in mathematical physics. 

An in depth critical analysis of the work of Milek and Seba was also under¬ 
taken; I highlighted a number of misconceptions and errors in the work. I have 
also highlighted that, even with Bourget’s demonstration of the existence of a con¬ 
tinuous spectral component of the Floquet operator for J > 3, Milek and Seba’s 
work still is not fully justified. Using Vinogradov’s conjecture, one can then jus¬ 
tify Milek and Seba’s result. A direct proof of Combescure’s conjecture or a proof 
of Vinogradov’s conjecture, allowing my work to bridge the gap, remains desir¬ 
able. 
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In this thesis I have presented new results on the elassifieation of the speetrum 
of the Floquet operator for a elass of kieked Hamiltonian systems. 

Before presenting the main results, I reviewed the fields of elassieal and quan¬ 
tum ehaos (Chapter 2) and presented a detailed investigation of the links between 
quantum dynamies and speetral analysis (Chapter 5). A point of ambiguity in the 
physies literature was also identified and eommented on. 

In Chapter 4 I provided a eonceptual introduetion to the mathematieal field of 
speetral and funetional analysis, highlighting the physieal meaning behind mueh 
of the basie mathematieal building bloeks. The aim was not only to lay the ground¬ 
work for the material in Chapter 6, but to make this important field of mathematies 
more aeeessible to physieists. An appreeiation of funetional and speetral analysis 
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Conclusions and further work 


provides one with a deeper understanding of basie quantum meehanies. 

In Chapter 6 1 developed a number of unitary equivalent theorems to those well 
known in the self-adjoint theory (e.g., the Putnam-Kato theorem) and successfully 
applied them to show that the spectrum of the Floquet operator remains pure point 
(given that U = has pure point spectrum) when the perturbation is suit¬ 

ably constrained. It should be stressed that the result is non-perturbative. I also 
obtained this result in a more straight forward, but less general, way by extending 
the work of Combescure. This was presented in the early parts of Chapter 7. 

In the remaining sections of Chapter 7 I extended the results on the emergence 
of a continuous spectrum of the Floquet operator to rank-A^ perturbations and 
investigated Combescure’s conjecture that the eigenvalue sequence for the unper¬ 
turbed Hamiltonian, Hq, should not affect the results in a significant way. Re¬ 
viewing the work of Bourget and linking it to a number-theoretic conjecture put 
forward by Vinogradov, I showed that if one could solve Vinogradov’s conjecture, 
then the essence of Combescure’s conjecture would follow. For any eigenvalue 
sequence for Hq described by a polynomial, a continuous component will emerge 
under realisable conditions for the rank-A^ operator perturbation. 

8.1 Future directions 

The clearest loose end in the work presented in this thesis would be to find a proof 
of Vinogradov’s conjecture on the estimation of Weyl sums. Without a proof, it 
seems reasonable to conclude that the work of Milek and Seba on kicked rotors 
and the existence of a continuous component cannot be properly justified. Having 
stood for over fifty years, a proof or disproof seems unlikely to turn up any time 
soon. 

When considering open questions of physical, rather than mathematical, in- 
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terest, a different direetion is seen for extending this work. It is broadly ac¬ 
knowledged that moving beyond the closed system dynamics governed by the 
Schrodinger equation is of great use in analysing chaos in quantum systems and 
investigating the quantum-classical correspondence. Thus, it would be interest¬ 
ing to allow for environmental interactions and to study the behaviour of kicked 
systems in so called “open quantum systems”. There is already a great deal of 
research effort in open systems and specifically, the analysis of kicked systems in 
such environments. Many of the results look promising. Applying the techniques 
developed in this work would be interesting. 

Returning to the closed system dynamics, an extension of the work presented 
in Chapter 6 to infinite N would be valuable. Similar work for the Hamiltonian 
system was done in the later part of Howland’s paper [66]. As > cx) a number 
of the tools and arguments used for finite N clearly collapse and one must be 
rather cautious. The physical implications of infinite N systems are also open to 
investigation. 
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